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We provide a framework for analyzing the problem of interacting electrons in a ballistic quantum 
dot with chaotic boundary conditions within an energy Et (the Thouless energy) of the Fermi 
energy. Within this window we show that the interactions can be characterized by Landau Fermi 
liquid parameters. When g, the dimensionless conductance of the dot, is large, we find that the 
disordered interacting problem can be solved in a saddle-point approximation which becomes exact 
as g — > oo (as in a large- N theory). The infinite g theory has two phases as a function of the Landau 
parameter u m in a channel with angular momentum m: A weak-coupling phase where constant 
charging and exchange interactions dominate the low-energy physics, as in previous "Universal 
Hamiltonian" treatments, and a strong-coupling phase characterized by the same order parameter 
as in the Pomeranchuk transition in clean systems (a spontaneous interaction-induced Fermi surface 
distortion), but smeared and pinned by disorder. Thus, both interactions and disorder are crucial 
to the existence of these phases. At finite g, the two phases and critical point evolve into three 
regimes in the u m — 1/g plane - weak- and strong-coupling regimes separated by crossover lines 
from a quantum-critical regime controlled by the quantum critical point. In this, the first of a two 
part series, we focus on the consequences of this picture for Coulomb Blockade experiments. We 
employ analytical and numerical methods to predict the statistics of single-particle levels, Coulomb 
Blockade peak spacings, conductance peak heights and quasiparticle widths. We show that in the 
strong-coupling and quantum-critical regions, the quasiparticle acquires a width of the same order 
as the level spacing A within a few A's of the Fermi energy due to coupling to collective excitations. 
In the strong coupling regime if m is odd, the dot will (if isolated) cross over from the orthogonal 
to unitary ensemble for an exponentially small external flux, or will (if strongly coupled to leads) 
break time-reversal symmetry spontaneously. For any m, the peak spacing distribution becomes 
broader than expected in previous works and even has support at negative values, which in turn is 
correlated with small peak heights. Ballistic/chaotic quantum dots afford us unrivalled theoretical 
and experimental control over the problem of simultaneous disorder and interactions due to the 1/g 
expansion and our ability to vary disorder and interaction much more readily than in the bulk. 
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I. INTRODUCTION 

The problem of treating electronic interactions in 
mesoscopic systems brings together two very interesting 
subfields of condensed matter physics. On the one hand, 
bulk systems with interactions and disorder can show un- 
expected phenomena 1 , and are the subject of ongoing 
and vigorous investigation. On the other hand, meso- 
scopic systems 2 can show behavior which is not present 
in bulk systems. Two examples are the oscillations as 
a function of gate voltage V g in the tunnelling conduc- 
tance of a quantum dot (QD) weakly coupled to the leads 
(the Coulomb Blockade (CB) regime 3-7 ), and persistent 
currents 8 in small metallic rings subject to a Aharanov- 
Bohm flux. 

In recent months we have developed a formalism 
that can tackle the problem of mesoscopic, disordered, 
strongly correlated systems 9,10 . We showed that there 
exists a controlled and complete solution to this problem 
provided the quantum dot (or ring) is ballistic, meaning 
that the disorder comes from chaotic collisions with the 



walls of the mesoscopic structure, and not from impuri- 
ties inside it. The solution becomes exact in the limit 
when the dimensionless conductance g of the quantum 
dot becomes large 9,10 . In this paper we focus on the de- 
tails of our approach in CB regime, while Part II 11 of 
this paper will analyze the effects on persistent currents 
(a short report of which has appeared 12 ). Experimen- 
tal Coulomb Blockade samples 5-7 are ballistic and do 
have fairly large conductances (g « 5 — 20). Our predic- 
tions are directly applicable to such samples even if the 
interactions between electrons are strong, which can be 
achieved by lowering the electron density below that of 
present samples. 

Let us first look at some important length, time, and 
energy scales to set up the problem. We will confine 
our analysis to two dimensions, since experimental sam- 
ples are made by confining a two-dimensional electron gas 
(2DEG) laterally by means of gates. We will focus atten- 
tion on ballistic QD's for which the bulk mean free path 
I is much larger than system size L. Ignoring the (weak) 
coupling to the leads produces sharp single-particle levels 
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in the QD whose mean spacing is A. A famous conjecture 
by Bohigas, Giannoni, and Schmidt 13 states that all lev- 
els within a Thouless energy Ep — Tivp/L of each other 
have correlations controlled by Random Matrix Theory 
(RMT) 14 . The analogous conjecture has been proved 
for the case of diffusive QD's 15 , and there are ongoing 
efforts to prove it for the ballistic case 16 . We will be in- 
terested in applying RMT to states that lie within Ep 
of the Fermi energy. Of particular interest to us are two 
RMT ensembles 14 , the Gaussian Orthogonal Ensemble 
(GOE) which applies to time-reversal invariant systems, 
and the Gaussian Unitary Ensemble (GUE) which applies 
to systems in which time-reversal is broken by an external 
magnetic flux. We will assume that spin-orbit coupling 
is negligible throughout this work. Also, because of a 
variety of factors, the coupling of the external magnetic 
field to the spin is small compared to the orbital effects 
in GaAs, the standard material used to make quantum 
dots. To simplify matters, we will set the spin-B-field 
coupling to zero in what follows. The dimensionless con- 
ductance g (also called the Thouless number) is defined 
by the ratio g = Ep/ A. Two energies inherited from the 
bulk 2DEG are the bandwidth Et and the Fermi energy 
Ep. The strength of interactions in the bulk 2DEG is 
characterized by the dimensionless number r s = a/ao, 
where a = 1/ y/irp is the typical distance between neigh- 
boring particles (p is the density) while ao = Ti 2 /me 2 is 
the Bohr radius. The band effective mass must be used 
in computing ao as must the dielectric constant. Since 
the system is finite, capacitive effects produce a charging 
energy Uq = e 2 /C, where C is the capacitance. 

The experiments that contributed greatly to our un- 
derstanding have been on the the zero-bias conductance 
of a QD weakly coupled to the leads as a function of V g 
at T = 0. At a generic value of V g the ground state 
has a definite number of particles N and energy £m- If 
the chemical potential fi = £n+i — £n — cteV g (a is a 
geometry-dependent "leverage" factor ) the free energies 
of the N and N + 1-particle states are degenerate, and a 
tunneling peak occurs at zero bias. Successive peaks are 
separated by the second difference of £n, called A2, the 
distribution of which is measured 5 ' 6 . 

To analyze the problem theoretically, it is simplest 
to consider a theory within the Thouless shell around 
E F (defined by \e - E F \ < E T /2). In this shell all the 
statistical properties of the single-particle energies and 
wavefunctions can be obtained from RMT. The generic 
Hamiltonian in this shell can be written as 

a a(3j8 

where e a are single-particle levels that obey RMT statis- 
tics, have a mean spacing A and range from —g A/2 to 
g A/2, and V a ^s represents a two-body interaction. In 
the following we will supress spin for simplicity, pointing 
out how its restoration modifies various results. Differ- 
ent forms of Vap-yg turn out to represent very different 
physics. 

The simplest model for interactions in a QD has a 
constant charging energy 3 ' 4 Uq: 



Hu = J2^cl s c a ,s + Y^ 2 ( 2 ) 

a.s 

which corresponds to the choice V a pp a = Uq = uq A, (in- 
dependent of a and /3 with all other couplings not of this 
form vanishing). This model predicts a bimodal distri- 
bution for A2: Adding an electron above a doubly-filled 
(spin-degenerate) level costs Uo + e, with e being the en- 
ergy to the next single-particle level. Adding it to a singly 
occupied level costs Uo- While the second contribution 
gives a Dirac delta-function peak at Uo, the first contri- 
bution is the distribution of nearest neighbor level sepa- 
ration e, which is known 14 from RMT to have a width of 
the order of the mean single-particle level spacing A. The 
puzzle that stimulated recent theoretical developments in 
the field is that numerics 1 ' 18 and experiments 5,6 produce 
distributions for A 2 which do not show any bimodality, 
and are much broader. 

One interesting option is to consider U Q( 3 7< 5 as indepen- 
dent gaussian variables 19 . Our work is more closely re- 
lated to another, the Universal Hamiltonian 20 ' 21 , wherein 
( for the spinful case) a term coupling to total S 2 and a 
Cooper coupling also appear: 

Hu= E z a clsc a . s + f N 2 - ±S 2 

a, s 

+ a(e4,t c L)(Ec/ 3 4C /3 , t ) (3) 

The rationale for the above choice of interactions is the 
following: As we will see below, only the "diagonal" ma- 
trix elements in which the indices are pairwise equal sur- 
vive the ensemble average and the variances of all the 
matrix elements are small 

A 2 

< V 2 0lS >-< V a0lS > 2 ~ -j- (4) 

where the indices now describe orbital and spin variables. 
Thus for large g, the interaction may be approximated 
by its ensemble average, which is just the Universal 
Hamiltonian. Apart from this motivation, calculations 
based on the Universal Hamiltonian 20 ' 21 have proven very 
successful at describing experiments 5 on quantum dots 
with small r s w 1 at a quantitative level, once exper- 
imental noise is subtracted out and finite-temperature 
effects 22-24 are taken into account. 

Our work deals with situations where r s is a lot larger 
than in current samples, in which case other unexpected 
phenomena appear possible. Our choice of V a fjjs fol- 
lows from the assumption of Landau Fermi liquid in- 
teractions at the Thouless energy 9 ' 10 , which we justify 
using renormalization group arguments. The Landau 
interactions 25 ' 26 are parametrized by couplings u m in ev- 
ery angular momentum channel m. For the spinful case 
the set of parameters is doubled, with one for the charge 
and one for the spin. Keeping only UQ har9e and u s pm 
yields the Universal Hamiltonian (without the Cooper 
term) . The retention of all the Landau parameters leads 
to the emergence of strong-coupling regimes in m 7^ 



2 



angular momentum channels 9,10 , which exhibit a sponta- 
neous deformation of the Fermi surface (suitably smeared 
out and pinned by disorder) as well as the possibility 
of time-reversal violation. The transition was originally 
discovered in the clean bulk limit by Pomeranchuk 27 . 
While these are true second-order quantum phase transi- 
tions in the limit g — > oo, they are replaced by sharp 
crossovers for finite but large g. This paper presents 
in detail the physics of these new phases, the transi- 
tions/crossovers between the weak- and strong-coupling 
regime via the quantum-critical regime and their atten- 
dant signatures. We concentrate on charge-channel in- 
stabilities in the spinful system; spin-channel instabilities 
necessarily involve the exchange coupling J and will be 
the subject of future work. The purpose of this set of two 
papers is to supply all the details left out of the previous 
brief presentations 9,10,12 , as well as to present numerical 
corroboration of the earlier analytical results, many new 
results, and experimental signatures. 

The plan of the paper is as follows. In section II 
we present the renormalization group (RG) argument 28 
for starting with Landau Fermi liquid interactions in the 
Thoulcss band, with all the assumptions and caveats. In 
particular, we will show why these assumptions are rea- 
sonable for ballistic QD's but might fail for QD's deep 
in the diffusive limit 29 . In Section III we present a fur- 
ther RG analysis inside the Thoulcss shell (modeled after 
Ref. 30 ), this time integrating out exact eigenstates of the 
disordered single-particle Hamiltonian. To leading order 
in 1 1 g we will see that the RG flow to one-loop order im- 
plies the phase transitions mentioned above. However, 
some questions about the effects of higher loops are left 
open, as is the very nature of the strong-coupling phase. 
These issues are resolved in Section IV, where we show 
that the problem can be mapped on to a large- N theory, 
with g playing the role of N. * The largeness of g allows 
us to fully control the calculation at at strong-coupling 
and not only confirm the phase transition discovered in 
perturbative RG 9 but also to probe the strong coupling 
phase in considerable detail 10 . In the next and longest 
section (Section V) we explore the physical consequences 
of the large-g solution: The nature of the order parame- 
ter, the physical signatures of the strong-coupling regime 
and the crossover quantum-critical regime 31 , and some 
curious properties of the weak-coupling regime. Numeri- 
cal methods which are capable of accessing all regions of 
the phase diagram (except the quantum critical regime) 
corroborate the analytical methods and also complement 
them, for example in the computation of level statistics 
and the CB peak-spacing statistics in the strong-coupling 
regime. 

We predict that if the system is in the strong-coupling 



"This large- N nature has two manifestations, one which is 
best expressed by saying that a certain class of diagrams dom- 
inate over others, and another in which a large number (here 
g 2 ) appears in front of the effective action, showing that the 
saddle-point approximation is a good one. 



regime of an even m channel it will display a broaden- 
ing of the CB peak spacing distribution, correlations be- 
tween small peak-heights and small peak-spacings, and 
a diamagnetic persistent current (which is not expected 
theoretically 32,33 in non-superconducting materials, but 
is seen experimentally 8 ). 

In the strong-coupling regime of an odd m channel, in 
addition to all the above effects (but with a paramag- 
netic persistent current), the system can be mapped on 
to the Caldeira-Leggett 34 model of a particle in a double- 
well potential subject to ohmic dissipation. The solution 
to this model 35,36 shows that the m odd system can be 
driven through a quantum phase transition which sponta- 
neously breaks time-reversal symmetry by changing the 
coupling of the quantum dot to the leads. (This really 
is a phase transition since the quantum coherence of the 
dot with the infinite reservoir renders a sharp transition 
possible). 

In this section we will also show that while an ef- 
fective single-particle description works in the weak- 
coupling regime, it fails in the quantum-critical and 
strong-coupling regimes even at low energies of order a 
few A. In these latter regions we will make a connection 
to ideas of "Fock-space localization/delocalization" 37,38 , 
which captures the crossover of the spectral width of the 
quasiparticle to a Breit-Wigncr form as the quasiparticlc 
energy is increased. This long section ends with its own 
summary. In Section VI our conclusions are presented, 
as are discussions of how our approach might be useful in 
other problems involving disorder and interactions. Ap- 
pendices A and B contain some additional details. 



II. EFFECTIVE INTERACTIONS AT THE 
THOULESS ENERGY 

In this section we will give a brief introduction to Fermi 
liquid theory 25,26 in the bulk, describe its instabilities 27 , 
and see how it must be modified to account for disorder. 
We will argue that a hybrid Hamiltonian, with a single- 
particle part encoding the chaotic nature of the single- 
particle states, and a part with Fermi liquid interactions, 
is the most natural starting point for analyzing ballistic 
mesoscopic structures in the Thouless band. 



A. The Clean Limit in the Bulk 

It is well-known that in a clean two-dimensional bulk 
system, no matter what interactions one starts with in 
the bare Hamiltonian, provided no superconducting or 
density wave instabilities intervene, the effective inter- 
actions are of the Landau Fermi liquid form 25,26,28 in a 
sufficiently thin shell near the Fermi surface defined by 
|k| — k F < k max . This corresponds to an energy cutoff 
A = vpk max . The entire Hamiltonian in such a thin shell 
can be written for the spinless case as 

H FL = Y> (k) c t(k) c (k) + ^ £ u{9-8>) 

k k,k',q 
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x :ct(k-q)c(k) C t(k' + q)c(k') : (5) 

where £o(k) is the renormalized quasiparticle energy 
measured from the Fermi surface, c(k), c' (k) are canon- 
ical anticommuting fermion operators, N n is the single- 
particle density of states, the :: sign stands for normal- 
ordering (subtracting the average in the ground state), 
9, 9' are the angles of the two-dimensional vectors k, k', 
and u(9 — 9') is the (dimensionless) Fermi liquid in- 
teraction function. It is understood that k, k' are 
to be summed only over states in the shell, and that 
|q| < k max . The main features of the Fermi liquid inter- 
action are that u{9 — 9') does not depend on the radial 
magnitudes of k, k' but only on the angles, and that the 
interaction shows mostly forward scattering. Rotational 
invariance has been used to express u as a function of the 
difference of the angles. One then Fourier decomposes u 
to obtain 

oo 

u{6- 9') = M + u m cosm(9- 9') (6) 

m— 1 

where the u m are known as Landau parameters. The 
Fermi liquid interaction function also determines the en- 
ergy to add a quasiparticle with momentum k in the pres- 
ence of a background of excited quasiparticles/quasiholes 
specified by Sn(k') (which is the deviation in occupation 
from the ground state) 

e(k)^ (k) + i^^-fi')Mk') (7) 

and the Landau total energy functional 

£({Sn}) = ]T eo(k)Mk) + ^E u ( fl - 9')6n(k)6n(k>) 

k kk' 

(8) 

Landau's original derivation 25 ' 26 rests on phase space 
and adiabatic continuity arguments and also predicts 
a decay rate for quasiparticles of order s 2 /Ef, which 
has been neglected above. The same result can also 
be derived by integrating out high energy states (single- 
particle states far from the Fermi surface) in a renormal- 
ization group (RG) approach, as was shown by one of 
us a decade ago 28 . In this approach, at a given stage 
of the RG one has a theory with a certain cutoff k max . 
One then integrates out a thin "shell" of momenta or 
width Sk max , thereby obtaining an effective theory with 
a cutoff k max — Sk max . One demands that the new the- 
ory give the same answer for all physical Green's func- 
tions in the low-energy sector as the theory with cutoff 
kmax- To achieve this the coupling constants must flow 
as one changes the cutoff. If a certain coupling constant 
increases as high-energy states are integrated out, it will 
end up dominating the low-energy physics no matter how 
small it was initially. Such a coupling is called "relevant" . 
In the opposite case a coupling may shrink as high-energy 
states are integrated out, in which case it is unimportant 
for low-energy physics, and is called "irrelevant" . Cou- 
plings which do not flow are called "marginal" . It turns 



out that all the Landau parameters arc marginal for a 
clean system in the bulk, whereas all other types of cou- 
plings, with one exception are irrelevant. The exception 
is the coupling in the BCS channel, i.e., between parti- 
cle of opposite momenta. These are irrelevant if repulsive 
and relevant if attractive. Since no superconductivity has 
ever been detected in the GaAs 2DEG's which are the 
basis of the ballistic quantum dots we focus on, we will 
assume that the BCS instability is absent. 

Thus, the Landau theory defines a fixed point (in fact 
a whole class of fixed points) for the clean electron gas. 

B. Fermi Liquid Parameters for the Spinful Case 

In the spinful case the Fermi liquid interaction function 
is to be thought of as a matrix in spin space 26 . Consider 
the energy to add a particle in a momentum state k with 
a density matrix p ss i (k) . This energy is 

e {p} (k) = Tr(e(k)p(k)) = ^ e ss ,p s , s (9) 

ss' 

The energy matrix e ss ' (k) depends on the occupations 
and spin states of other quasiparticles, or in other words, 
the density matrices dn SlS > i (k') in the following way 

£ S s'(k) = £o,ss'( k ) + 

^ E u ss ,, Sls[ (6-6')5n Sls[ (V) (10) 

si s^ ,kk' 

with the corresponding total energy functional 
£({Sn}) = E £o, ss '(k)5n s / s (k) 

ss'k 

+ 2^ E u ss ,, sis[ (e-6')dn ss ,(k)5n sis[ (k') (11) 

ss' s-is^ ,kk' 

We will restrict ourselves to the case with spin- 
rotation-invariance, which leads to a restricted form for 
the interaction function 26 

u«',. ia i (0 - 0') = $(9 - 9')S SS ,S S1S[ + Z(9 - 9')f ss , ■ f SlS , 

(12) 

where the f = (t^ x \t^ v \t^) are the Pauli spin ma- 
trices, and the $ and Z are the Fermi liquid interaction 
functions in the charge and spin channels respectively. 
As usual, these can be Fourier expanded to obtain the 
charge and spin channel Landau parameters 

oo 

$(9 -6') = $ + E *m cos m(9 - 9') (13) 

ra— 1 

oo 

Z{9 -6') = Z + E z m cosm(9 - 9') (14) 

m— 1 

The $ m and Z m are identical to the parameters used 
by Pines and Nozieres 41 , 

F m = <& m (15) 

F m = Zm (16) 
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The Landau parameters depend on the relative 
strength of the interaction and kinetic energies, com- 
monly characterized by the dimensionless number r s — 
a/ do, where a = is the typical distance be- 

tween neighboring particles (p is the density) while a = 
h 2 /me 2 is the Bohr radius. The band effective mass must 
be used in computing ao as must the dielectric constant. 
The dependence of the Landau parameters on r s has been 
the subject of investigation by quantum Monte Carlo 
methods 42 . For the largest value investigated r s = 5, 
the results in our normalization (differing from that of 
ref. 42 by a factor of two) are 



$0 = 


- 1.85 


(17) 


Za = 


- 0.25 


(18) 


*i = 


0.06 


(19) 


Zi = 


- 0.135 


(20) 


$ 2 = 


- 0.25 


(21) 


z 2 = 


0.16 


(22) 



in the spin channel for illustration. Now Sn must be 
a density matrix in spin space, and we can choose, for 
example, 

Working through the energy functional of Eq. (11) in 
the same way as before one obtains for m ^ 

where the extra factor of two inside the brackets com- 
pared to Eq. (25) comes from the extra trace over the 
two-dimensional spin space, combined with the fact that 
in our normalization, Nq is the spinless single-particle 
density of states. Thus, in the spinful case, with our 
normalizations, the instability occurs for <& m , Z m < — 1. 
From the values quoted 42 in the last section, we see that 
the electron gas in the clean bulk limit is quite far from 
these instabilities even for t s — 5. 



C. Instabilities of the Clean Fermi Liquid 

It is well-known that the Fermi liquid is unstable to- 
wards the introduction of an attractive coupling in the 
Cooper channel which leads to a gapped superconducting 
ground state. The Fermi liquid also has other instabil- 
ities for certain values of the Fermi liquid parameters, 
a fact first pointed out by Pomeranchuk 27 . There has 
been a revival of interest in the Pomeranchuk transition 
in clean bulk 2D systems recently 43 ' 44 . 

Let us first consider the spinless case. We consider a 
deformation of the Fermi surface by an amount r{9) in 
the direction 9. To be precise, we consider the following 
deformation 



Sn(k) = 



1, < e (k) < r(0), r > 
-1, r{9) < e (k) < 0, r < 



(23) 



Using the replacement J2k ^> No J de we can calculate 
the energy cost for this deformation using the Landau 
energy functional, Eq. (8), 



iVo 



2tt 271 

f d9 r{9) 2 1 r 
J 2^^~ + 2j 



2ir 

+ 1 J U^-^'MW) (24) 
o 



To isolate a particular Fermi liquid channel of angular 
momentum m^Owe choose r(9) — r$ cos m9 to obtain 
the total energy for m ^ 



^ = i ( i 

N 4 V 



./2) 



(25) 



It is clear that for u m < —2 the energy is negative, 
indicating an instability of the undeformed Fermi liquid. 
This is the Pomeranchuk instability 27 . 

In the spinful case the "deformation" can be in the 
charge or the spin channel. Let us consider the instability 



D. The Disordered Bulk 

Let us now examine the changes that occur when 
elastic impurity scattering is taken into account in the 
bulk. The RG method can be conceptually generalized 
for this case as well. An important energy scale here is 
E T = K/t = Tivf/1, where I is the mean free path due 
to impurities and r the mean free time. In a time-scale 
r momentum states get scattered by impurities imply- 
ing that E T is the spectral broadening of a momentum 
state. Since impurity scattering is elastic, an exact eigen- 
state of the disordered single-particle Hamiltonian with 
energy e can be expressed as a superposition of momen- 
tum states of energy e with a spread of E T . Therefore for 
high energies £> £ T integrating out a momentum states 
of thickness E T is roughly the same as integrating out a 
corresponding shell of exact energy eigenstates. This is 
equivalent to the statement that states of high energy 
have negligible correlations with low-energy states. 

One can now ask how the interaction coupling con- 
stants flow in the disordered case. Consider first the case 
of a microscopic short-ranged interaction of range £, such 
as in He 3 . The time taken by a single fermion-fermion 
collision is then roughly £/vf- As long as this time is 
much smaller than r, the collisions happen far from an 
impurity, and therefore conserve total momentum. Thus, 
despite the fact that momentum is not a good quantum 
number for single-particle states, interactions continue to 
conserve momentum. Since integrating out momentum 
states is roughly the same as integrating out exact energy 
eigenstates, the conclusion is that up to an energy cutoff 
A of the order of magnitude of E T the RG flow should 
proceed much as in the clean system, resulting in Fermi 
liquid interactions. 

Once the cutoff reaches E T integrating out momen- 
tum states is no longer approximately equivalent to in- 
tegrating out exact eigenstates of the disordered single- 
particle Hamiltonian, since now the wavefunctions of the 



5 



states being integrated out have non-negligible correla- 
tions (induced by disorder) with states being kept. In the 
disorder-averaged version of the theory 45 ' 46 , these cor- 
relations manifest themselves as collective modes such 
as diffusons and Cooperons with singular propagators. 
In the diffusive limit, these propagators drive the RG 
flow 45 ' 46 , resulting in a runaway flow of the s-wave triplet 
coupling Z towards strong coupling. The resulting state 
seems to be weakly ferromagnetic 47 , but is currently not 
fully understood. 

Now consider a model with Coulomb interactions. Up 
to an energy scale of E T a good model for the effective 
interactions is the Thomas-Fermi static screened interac- 
tion vpp(q) = 2ire 2 /(q + qpp), with a screening wavevec- 
tor qpp ~ flo 1 m t w0 dimensions (do is the effective Bohr 
radius in the material) . By following the previous logic, it 
is clear that as long as 1/qppl is small, electron-electron 
collisions will conserve total momentum, leading to Fermi 
liquid interactions near Ep. However, at very low en- 
ergies, deep in the diffusive regime {lot -C 1, ql -C 1) 
this interaction can get unscreened 29 . The resulting low- 
energy unscreened interaction leads to a breakdown of 
Fermi liquid theory in two dimensions at the most fun- 
damental level, namely, the quasiparticle lifetime broad- 
ening becomes comparable to its energy 29 . It is also clear 
that in this regime, the electron-electron collision time is 
long compared to r, and therefore momentum will not 
be conserved in an electron-electron collision. 



E. Ballistic Quantum Dots 

Now let us examine the case of ballistic QD's, which 
is conceptually the same as considering electrons in a 
cavity with hard walls. Here since the role of impurity 
scattering is being played by scattering from the walls 
of the cavity, L plays the role of /. Therefore the Thou- 
less scale Ep is roughly the same as E T . Since we know 
that the interaction is screened at this energy scale, we 
can conclude that Fermi liquid interactions are valid at 
this scale. Can Fermi liquid interactions become invalid 
at lower energies due to Altshuler-Aronov-like effects 29 ? 
The answer is no: Due to the finite size of the cavity, 
all momenta can only be defined up to an uncertainly of 
2-k/L. Since q can only be defined up to order l/L, one 
can never achieve the condition qL <C 1. The deep diffu- 
sive regime where interactions get unscreened is inacces- 
sible in ballistic QD's, and is superceded by the Random 
Matrix (RMT) regime which we will analyze in detail in 
the next section. The conclusion of this set of arguments 
is that Fermi liquid interactions are natural for ballis- 
tic QD's, while they may not be for strongly diffusive 
{I < L) QD's. 

Let us now explicitly consider what the low-energy ef- 
fective Hamiltonian looks like for a ballistic QD in the 
Thoulcss band. We will first write it down, and then 
discuss it. 

H = ^2s a cl t c a + ^-^2u(e-9'):n k n^ : (28) 

a k,k' 



where a are labels for the exact eigenstates {g in 
number) of the singlc-quasiparticlc Hamiltonian encod- 
ing the chaotic boundary scattering, A (the mean level 
spacing) is the inverse density of spinless states 1/iVo, 
nk = c^(k)c(k), and it is understood that the sum over 
k, k' goes over g values near the Fermi surface. (We are 
using nk instead of Sn^, the difference being absorbed in 
a shift in the chemical potential.) The Hamiltonian is 
expressed in a hybrid basis, with the connection between 
the a and the k basis being given by the wavefunctions 
of the exact eigenstates a in the momentum basis 4> a (k) . 
The statistics of the energies e a and wavefunctions <\> a (k) 
are assumed to be controlled by RMT. Note that the mo- 
mentum k is uncertain up to 2ir/L, and it really repre- 
sents a patch on the Fermi surface. To be more precise, 
we can define the state we label by k in the dot as 

|k)=^| a )(a| e ^)=^|a)^(k) (29) 

a a 

For convenience we will choose the k labels to be equally 
spaced (separated by an angle 2n/kpL = 2n/g) on the 
Fermi circle. This set of states clearly satisfies all the 
boundary conditions of the dot, and is unitarily related to 
the set of disorder eigenstates \a). In the original Fermi 
liquid Hamiltonian (Eq. (5)) there was an additional sum 
over q. Here, since q is also uncertain up to 2ir/L and 
only scattering within the Thoulcss shell is allowed, the 
sum over q is eliminated. Since the k take on g discrete 
values, Landau parameters with angular momentum m 
of order g are not sensible in this theory. Our focus will 
be on small m of order 1, because these are the channels 
in which the Landau parameters are expected to have 
large magnitudes. 

An interesting fact about ballistic/chaotic quantum 
dots is that the ratio Ep/Ep is related to g: 

E F h 2 k 2 F /2m* k F L 
Ep-^P-Jl-- — ^ 9 (30) 

Recall that Fermi liquid theory becomes applicable only 
in an energy shell of size much smaller than Ep around 
the Fermi surface. Therefore the largeness of g is also 
a sufficient condition for the validity of the Fermi liquid 
interactions. 

To summarize, in writing down the Hamiltonian of Eq. 
(28), we are appealing to the fact that whereas single par- 
ticle momentum is not conserved due to collision with the 
walls, the momentum of a pair is conserved during a colli- 
sion if it takes place over a time scale shorter than L/vp, 
the time to bounce off the wall. The interacting part of 
the Hamiltonian continues to have the same form as in 
the clean limit even in the presence of chaotic boundary 
scattering. 



F. Recovering the Universal Hamiltonian 

Let us see how the Universal hamiltonian of Eqn. (3) 
emerges, starting with the spinless case. One can use the 
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(j> a (k) to express the entire Hamiltonian in the a basis in 
the form of Eq. (1), with 



kk' 



'*(k)^(k')-0*(k')^(k) 

x [0 7 (k')0 5 (k)-0 7 (k)^(k')] 



(31) 



These coefficients V a p^s will vary from sample to sam- 
ple. Let us focus on terms that survive the ensemble av- 
erage, which here reduces to a Random Matrix Theory 14 
average. Let us concentrate on the GOE, where all the 
wavefunctions in the real-space basis 4> a (r) can be made 
real. Since we are working in an approximate momen- 
tum basis labeled by the g values of k, the reality of the 
real-space wavefunctions translates to 



6 Q (k))* =0 a (-k) 



(32) 



The fundamental RMT wavefunction correlator is 

^a/3<5kk' 



< C(k)Mk') >= 



fJ 



(33) 



To find the average of V a 0^s we need the four-point 
correlator. Apart from exceptional cases when all the 
subscripts and/or all the arguments of the wavefunctions 
are equal, the four-point wavefunction correlator in the 
GOE can be written to leading order in 1/g as 



< C(ki)^(k 2 )0 7 (k 3 )0 5 (k 4 ) > 



9 2 
9 2 

ot &y S <5k! , - k 2 <5k 3 , - k 4 



(34) 



It is seen that only matrix elements in which the in- 
dices a/3j5 are pairwise equal survive disorder-averaging, 
and also that the average has no dependence on the en- 
ergy of a(3jS. In the spinless case, the first two terms on 
the right hand side make equal contributions and produce 
the constant charging energy in the Universal Hamilto- 
nian of Eq. (3), while in the spinful case (to be discussed 
at length in Section IIIB) they produce the charging and 
exchange terms. The final term of Eq. (34) produces the 
Cooper interaction of Eq. (3). 

Finally, one can explicitly calculate the variances 



{V a p 1 s) 



A* 



rn — 1 



(35) 



and see that they are small. Note that uo does not con- 
tribute to the fluctuations between different disorder re- 
alizations. 



III. INSTABILITIES OF THE UNIVERSAL 
HAMILTONIAN: RENORMALIZATION GROUP 
TREATMENT 

The Universal Hamiltonian 20 ' 21 is obtained by replac- 
ing the interaction by its ensemble average in the full 
Hamiltonian, arguing that sample-to-sample fluctuations 




FIG. 1. Feynman diagrams for the full four-point ampli- 
tude r a/ 3js- 



in the interaction matrix elements are small. At this 
point one can ask when small terms in the Hamiltonian 
can be safely discarded, and when they cannot 9 . As long 
as one is interested in low-energy properties, the RG is 
the perfect tool to answer this question. One simply in- 
tegrates out high energy states, and looks at the fate of 
the originally small couplings. If they are irrelevant in the 
RG sense, discarding them is justified, while if they are 
relevant they dominate the low-energy physics regardless 
of how small they were initially. 

Two of us carried out just such an RG analysis 9 , to 
the presentation of which this section is devoted. We 
will first consider the spinless case, and then make some 
remarks about the spinful case. We will close this section 
with a list of questions left unanswered by the RG. 



A. RG for the Spinless Case 

We start with the Hamiltonian for spinless fermions 
in a ballistic QD, Eq. (28), which we reproduce here for 
convenience 



n k n k / 



(36) 



k,k' 



The strategy for carrying out RG in a finite system 30 
is the following: (i) Since we are in the Thouless shell we 
cannot integrate out momenta, but must integrate out 
exact eigenstates of the chaotic single-particle Hamilto- 
nian. (ii) After integrating out some exact eigenstates 
at a given stage in the RG we have g' = ge~^ levels left 
(here £ is called the flow parameter of the RG). At this 
stage we compute a scattering amplitude r Q ^ 7< 5 for the 
process in which two fermions originally in states a(3 are 
scattered into states 7<5. This scattering can proceed di- 
rectly through the vertex V r Q , i a 7 i(^), or via intermediate 
virtual states higher order in the interactions, which can 
be classified by a set of Feynman diagrams, as shown in 
Figure 1. All the states in the diagrams belong to the g' 
states kept, (iii) We demand that the entire amplitude be 
independent of g' , meaning that the physical amplitudes 
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should be the same in the effective theory as in the orig- 
inal theory. This will lead to a set of flow equations for 
the V a p~/S- In principle this flow equation will involve all 
powers of V but we will keep only quadratic terms (the 
one-loop approximation). Then the diagrams are limited 
to the ones shown in Figure 1, leading to the following 
contributions to the scattering amplitude T al 3 7 s 

,N F (v)-N F (f>) 



+ £' 



£ 



VappvVv^S (37) 



, l-N F (p)-N F {u) 



where the prime on the sum reminds us that only the g' 
remaining states are to be kept and Np(a) is the Fermi 
occupation of the state a. We will confine ourselves to 
zero temperature where this number can only be zero or 
one, but the extension to finite temperature is straight- 
forward. Also recall that the matrix element Vap^s now 
explicitly depends on the RG flow parameter £. 

Now we demand that upon integrating the two states 
at ±g'A/2 we recover the same Ta^s- Clearly, since 



d 



d£ 9 Sg' 



(38) 



The effect of this differentiation on the loop diagrams 
is to fix one of the internal lines of the loop to be at the 
cutoff ±g'A/2, while the other one ranges over all smaller 
values of energy. In the particle-hole diagram, since \x 
or v can be at +g'A/2 or —g'A/2, and the resulting 
summations are the same in all four cases, we take a 
single contribution and multiply by a factor of 4. The 
same reasoning applies to the Cooper diagram. Let us 
define the energy cutoff A = g'A/2 to make the notation 
simpler. Since we are integrating out two states we have 
Sg' = 2 



= 



dV 



M E 



,N F (v)~N F (p) 



function, Eq. (31). Since there is a product of two V's 
in each loop diagram, and each V contains 4 terms, it is 
clear that each loop contribution has 16 terms. However, 
it turns out that only certain terms contribute to lead- 
ing order in the large-g limit, while others do not. We 
will focus on one of each type to illustrate the difference. 
Let us first consider a term of type I in the particle-hole 
diagram, which survives in the large-g limit. Putting in 
the full wavefunction dependences (and ignoring factors 
other than g, g') we have 

g' A 2 £ 



+ 



%4 E 



fi—AyU 



, l-N F (ti)-N F (iy) 



The changed sign in front of the 1-loop diagrams reflects 
the sign of Eq. (38) 

So far we have not made any assumptions about the 
form of Vafjjs, and the formulation applies to any finite 
system. In a generic system such as an atom 30 , the ma- 
trix elements depend very strongly on the state being 
integrated over, and the flow must be followed numeri- 
cally for each different set af3^5 kept in the low-energy 
subspace. The crucial simplification in QD's comes from 
the fact that these matrix elements are random and 
controlled by RMT, and that the RG equation is self- 
averaging, as will be seen below. 

Let us go back to the properly antisymmetrized matrix 
clement defined in terms of the Fermi liquid interation 



A+|£„ 

£ £ «(0k - P )u(0p' - 0kO 

kk' pp' 

&(k)^(k')0 7 (k')&(k) 
0*(p)«(p')Mp)Mp') 



(40) 



The internal sum over v is self-averaging. While the 
most convincing way to show this is to compute its vari- 
ance, and see that it is of order 1 / ^fg times its average, 
this fact can be motivated in the following way: There 
is a sum over g' ^> 1 values of v with a slowly varying 
energy denominator, which makes the sum over v similar 
to a spectral average, which in RMT is the same as an av- 
erage over the disorder ensemble. A more sophisticated 
argument is presented in Appendix A. 

In RMT the wavefunction averages do not depend on 
the energy separations of the states, so the wavefunction 
average can be carried out separately. This average in 
the GOE for fj, ^ v and generic momentum labels is 



< <Au(Pl)4C(P2)<MP3) 



SiaS; 



14«23 



oi3<>: 



13"24 



V(P 4 ) > = 
51,-2^3,-4 



(41) 



The 1/g 2 term is the "naive Wick contraction" of 
leading-order RMT, but the 1/g 3 terms are necessary to 
maintain orthogonality between (i and v. The final term 
would be missing in the GUE. Substituting the correct 
momentum labels for the particle-hole diagram we see 
that the wavefunction average is 



(39) 



°pp' 

g 2 



l 



g° 



(42) 



Using the self-averaging shown in Appendix A, the 
first term of Eq.(42) forces p = p' in Eq. (40). For large 
g, using 



(43) 



we obtain a convolution of the two Fermi liquid functions 

oc 

5>(# k - 9 p )u(9 p -6v)= g(u 2 + u 2 m cosm(8 - 9')) 



m—1 



(44) 
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where we have reverted to the notation 9 = 9-^, 9' = 9^. 
In the second term of Eq. (42), the t> P .- P ' turns out to 
be subleading, while the other allows independent sums 
over p, p'. This means that only u contributes due to 
this term, which produces 



5>(0k-0 P M0p' -9k>) =9 2 ul 



(45) 



We still need to do the energy sum, which we replace 
by an integral 



n A 
^ A ATM ~ J ~AA + e 



e„=-A 



In 2 



(46) 



Feeding this into full expression for this contribution to 
the particle-hole diagram, we find it to be 



^-Aln2£ J2 u 2 m cosm(9 - 9') 

kk' \m=l 
^(k)^(k')0 7 (k')0 5 (k) 



(47) 



Notice that the result is still of the Fermi liquid form. 
It is also seen that the contribution from uq cancels, 
about which we will say more below. 

Now let us consider a contribution which is subleading 
in 1/g in the large-g limit. 

o 

- a' A 2 V i 

v— — A 

E E u{e k - e p )u(e p , - 0k') 

kk' pp' 

<£(p)^(k')^(k')&(k) 
^(k)C( P ')^(p)^(p') 



(48) 



Note that the momentum labels of </>* and (/>* have 
been exchanged compared to Eq. (40) and there is a 
minus sign, both coming from the antisymmetrization of 
Eq. (31). Once again we ensemble average the internal 
(i, v sum, the wavefunction part of which gives 



< 



;(k)C(p')^(p)^( P ') >= 



<5kp"5pp' _ ^kp + Jk.-p'flp.-p' 

g 2 g 3 



(49) 



It is clear that there is an extra momentum restriction 
in each term compared to Eq. (42), which means that 
one can no longer sum freely over p to get the factor of 
g in Eq. (44), or the factor of g 2 in Eq. (45). Therefore 
this contribution will be down by 1/g compared to that of 
Eq. (40). The general rule is that whenever a momentum 
label corresponding to an internal line in the diagram 
(here \i and v) is forced to become equal to a momentum 
label corresponding to an external disorder label (here 
a, (3, 7, or 5), the diagram is down by 1/g, exactly as in 
the 1/N expansion. 

Turning now to the Cooper diagrams, the internal lines 
are forced to have the same momentum labels as the ex- 
ternal lines by the Fermi liquid vertex, therefore they do 
not make any leading contributions. 



Finally, one can collect all the terms, and conclude 
that for m ^ 



du r< 



*(ln2)«* 



(50) 



It must be emphasized that uq does not flow. The 
physical reason for this is that it commutes with the one- 
body "kinetic" part, and therefore does not suffer quan- 
tum fluctuations. The above equation can be written in 
a more physically transparent form by using a rescaled 
(for m/0 only) variable 



e ^Un 



in terms of which the flow equation becomes 



du r 



(ln2){4 = /3(u m ) 



(51) 



(52) 



where the last is a definition of the /3-function. To un- 
derstand the meaning of u m we go back to Eq. (35), 
which expresses the variances of Va^s in terms of u m 's. 
We start with an ensemble of QD's with dimcnsionlcss 
conductance g. Suppose we wish to define a new ensem- 
ble with dimensionless conductance g' and Fermi liquid 
parameters u' m , and demand that all the statistical prop- 
erties of the matrix elements, including the variances of 
Eq. (35), be the same, we are led to 



u 



,'2 



(53) 
(54) 



It is clear that the u' m are exactly the u m . From 
Eq. (52) it can be seen that positive initial values of u m 
(which are equal to initial values of u m inherited from the 
bulk) are driven to the fixed point at u m = 0, as are neg- 
ative initial values as long as u m {^ — 0) > u* m = —1/ In 2. 
Thus, the Fermi liquid parameters are irrelevant for this 
range of starting values. Recall that setting all u m = for 
m/0 results in the Universal Hamiltonian 20 . Thus, the 
range u m > u* m is the basin of attraction of the Universal 
Hamiltonian. On the other hand, for u m (£ = 0) < u* m 
results in a runaway flow towards large negative values of 
u m , signalling a phase transition to a phase not pertur- 
batively connected to the Universal Hamiltonian. Recall 
that the Pomeranchuk instability 27 of the clean spinlcss 
Fermi liquid happens for u m < —2. The flow towards 
large negative values of u m suggests that the system is 
undergoing some kind of Pomeranchuk instability. This 
notion turns out to be correct, and will be made pre- 
cise in the next section. Note that there is a window 
— 2 < u m < u* n in which the clean bulk system is stable 
while the system in a ballistic QD is unstable. 

Let us summarize the results of the RG: (i) The RG 
equation is self-averaging for large g. (ii) The Fermi liq- 
uid form of the interactions is left unchanged to leading 
order by the RG flow, (iii) Each u m flows separately, 
while uq does not flow, (iv) The flow equations show an 
instability for u m (£ = 0) < u* m = —1/ In 2 in every Fermi 
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liquid channel, presumably towards a Pomeranchuk-likc 
state, (v) The critical point shows a "correlation length" 
exponent v = 1, calculated from the slope of the (3- 
function at the critical point. 



dfi 
du m 



= 1 



(55) 



is characterized by antisymmetry in spin, and therefore 
symmetry in the orbital labels of the incoming (and sep- 
arately, the outgoing) particles. The triplet channel is 
exactly the reverse, (ii) If s = s' then the interaction is 
forced to be in the triplet channel. However, if s ^ s' the 
interaction is a superposition of the singlet and triplet 
channels. 

On the above basis we identify the singlet and triplet 
channel interactions as 



B. Some Remarks on the Spinful Case 

The case with spin is similar conceptually but some- 
what more complicated technically than the spinless case. 
For a spin-rotationally-invariant Hamiltonian, the two- 
body interaction can be decomposed 25 ' 26 into a spin- 
singlet interaction function (9 — 6') and a spin-triplet 
interaction function u^(8 — 9'). The projectors for the 
two two-body spin states are 



1-T-f' 
4 

3 + r-r' 



(56) 
(57) 



where r and f 1 represent the Pauli spin operators for the 
two electrons (the numbers 1 and 3 in the numerators 
should be thought of as unit matrices in spin space). 
Writing the total interaction and comparing it to the 
Fermi liquid form (Eq. (12)) we find 



U (»)7?W + u (*)t?(*) = $ + Zt ■ f 



(58) 



where all Fermi liquid labels and spin-matrix indices 
have been suppressed. This leads to the following rela- 
tions 



t W + 3 U (*) 



,.(*) 



(59) 
(60) 



In order to carry out the RG, we define the interaction 
Hamiltonian as 



H 



V' 



2 v a/3yS'-as C f)s' C l s ' CSs 



(61) 



The one-loop scattering amplitude now reads as fol- 
lows 



af3jS 



V. 



+ E 



yss ys s 



V ss ,vs s 



- V s >s ' V ss ' — V V SSl V s Sl 
ftvfiy afivS ^s\ afivS ftvfij 

E l-N F (n)~N F (v) yss' yss' 



(62) 



In order to separate g into its singlet and triplet 
components, we need two facts: (i) The singlet channel 



y(s) 



y(t) _ 1 



v r 



afjfS ' {3a-y8 



u 



+ v r 



tl 



v, 



ri 



yU . 



yU 



yU , t/U 



(63) 
(64) 



and the inverse relations 



yU 



T/TT . 

V af3f8 ' 
' v a/3-yS 



y(t) 



V 



(t) 



(65) 
(66) 



In order to carry out the wavefunction averages one 
expresses these matrix elements as 



V 



vS, s = fE u^\e e') 0* (k)^(k') - 0* (k')^(k) 



±J2u {s) (9-9') 



kk' 



(k)^(k') + <&(k')^(k) 

x [0 7 (k')^(k)+0 7 (k)^(k')] 



(67) 



kk' 



x [0 7 (k')^(k)-0 7 (k)0i(k')] 



(68) 



Now it is straightforward to find the RG flows of the 
singlet and triplet Landau interaction functions. A non- 
trivial check is provided by the fact that there are two 
expressions for the triplet interaction, and their flows 
should be identical. These flows can be written in terms 
of the rescaled Landau parameters $ m = e^$ m and 
Z m = e~^Z m , and one obtains to leading order in 1/g 



d$ r 



d£ 

dZ m 



= -$ m -(21n2)*£, 
= -Z m -(2\n2)Zl 



(69) 
(70) 



Clearly, the charge and spin channels flow indepen- 
dently, and there is an instability at —1/(2 In 2) = 
—0.7213 ... in each channel. Recall that in the clean 
Fermi liquid the Pomcranchuk instability occured at —1 
in our normalization. Thus, there is a substantial win- 
dow where the clean system in the bulk is stable, but the 
ballistic system is unstable. 



C. Unanswered Questions in the RG 

While the RG makes it plausible that there is indeed 
a transition out of the phase controlled by the Univer- 
sal Hamiltonian, it leaves two important questions unan- 
swered. Recall that we have used a one-loop calculation 
of the scattering amplitude to define the RG flow. This 
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is usually sufficient when the dimensionless couplings u m 
are much smaller than unity. However, the critical point 
occurs when u m is of order unity. Can we trust the one- 
loop calculation at this point? In other words, suppose 
there were higher-order diagrams in the scattering am- 
plitude to leading order in 1/g. This would have gener- 
ated cubic and higher-order terms in the RG flow. For 
a sufficiently large value of the cubic term, the critical 
point identified in the previous subsections can disap- 
pear. In that case it would have been an artifact based 
on an invalid approximation. Secondly, even if we accept 
the phase transition as genuine, the RG gives us no way 
of calculating physical quantities in the strong-coupling 
phase. Both of these questions are answered in the next 
section, where we show that the one-loop flow is exact 
in the large-g limit, and that there is a controlled way of 
calculating all the properties of the strong-coupling phase 
in a large- N approximation. 



IV. LARGE- iV THEORY OF THE MESOSCOPIC 
POMERANCHUK REGIMES 

As the name suggests, large-iV theories typically in- 
volve interactions between N species of objects. The 
largeness of N renders fluctuations (thermal or quantum) 
small, and enables one to make approximations which are 
not perturbative in the coupling constant, but are con- 
trolled by the additional small parameter 1/N. Large- N 
approximations have found many uses in condensed mat- 
ter physics, but usually the value of N which character- 
izes the real model is quite small (2 or 3). We will find 
below that ballistic mesoscopic structures offer a realiza- 
tion of large- N theory where N — g. The g momentum 
states act as g "species" of fermions. Thus, one can make 
N as large as one wants modulo technological challenges. 
In fact, standard Fermi Liquid theory in a clean bulk 
system can also be formulated in this way 28 . 

There are two primary manifestations of the large- N 
nature of a theory: (i) The four-point scattering ampli- 
tude is dominated by iterated particle-hole diagrams, (ii) 
The effective action is dominated by its saddle-point. 

In this section we will demonstrate each of these man- 
ifestations of large-TV theory, and end by solving the 
strong-coupling theory in a saddle-point approximation. 



A. Diagrammatic Large-iV and the Exactness of 
One-Loop RG in the Large-g Limit 

To put the issue in context, let us consider the Gross- 
Neveu model 48 which is one of the simplest fermionic 
large-iV theories. This theory has N identical massless 
rclativistic fermions interacting by an attractive short- 
range interaction. The Lagrangian density is 




FIG. 2. The first few Feynman diagrams in the four-point 
scattering amplitude for the Gross-Neveu theory. 

Note that the propagator conserves the internal index, 
as does the interaction term. 

Figure 2 shows the first few diagrams in the expression 
for the four-point scattering amplitude (for particles of 
species i and j to scatter) in the Gross-Neveu theory. The 
"bare" vertex comes with a factor X/N. The one-loop di- 
agrams all share a factor A 2 /N 2 from the two vertices. 
The first one-loop diagram has a free internal summa- 
tion over the index k over N values, with the contribu- 
tion being identical for each value of k. Thus, this one- 
loop diagram acquires a compensating factor of N which 
makes its contribution of order A 2 /N, the same order in 
1/N as the bare vertex. However, the other one-loop di- 
agrams have no such free internal summation and their 
contribution is indeed of order 1/N 2 . Therefore, to lead- 
ing order in 1/N, one should keep only diagrams which 
have a free internal summation for every vertex, that is, 
iterates of the leading one-loop diagram, which are called 
bubble graphs. For later use remember that in the dia- 
grams that survive (do not survive), the indices i and j 
of the incoming particles do not (do) enter the loops. Let 
us assume that the momentum integral up to the cutoff 
A for one bubble gives a factor H(A,q ext ), where q ex t is 
the external momentum transfer at which the scattering 
amplitude is evaluated. To leading order in large- TV the 
full expression for the scattering amplitude is 



1 



A 



TVl + An(A, fet ) 



(72) 



Once one has the full expression for the scattering am- 
plitude (to leading order in 1/N) one can ask for the RG 
flow of the "bare" vertex as the cutoff is reduced by de- 
manding that the physical scattering amplitude T remain 
insensitive to the cutoff. One then finds 



dTij(q ext ) 



(IX 



— = X 



, dU{K,q ex t) 



(73) 



(71) 



where £ = ln(Ao/A) is the flow parameter of the RG. 
This equation shows that the one-loop RG flow is the 
exact answer to leading order in a large-N theory. All 
higher-order corrections must therefore be subleading in 
1/N. 
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Let us now turn to our theory, in which the role of the 
internal index is played by the momentum patch label k, 
which can take on g values. The role of the coupling con- 
stant A in our theory is played by Etu(8 — 9'). From the 
definition of the matrix elements in Eq. (31) we see that 
interactions occur in the Hamiltonian as Etu(6 — 9')/g 
in analogy with the X/N in the Gross-Neveu Lagrangian. 
Let us remark in this context that the most natural way 
to take the large-g limit is to keep Er fixed and decrease 
A. 

The major departure from the Gross-Neveu theory is 
that in our theory the propagator does not conserve the 
label k. The contributions of the leading and subleading 
one-loop diagrams were presented in Eqs. (40,48) respec- 
tively. The crucial fact which allows us to show the large- 
TV nature of this theory is the self-averaging of the inter- 
nal summations, which follows from the arguments of the 
previous section (and Appendix A). This self-averaging 
means that the internal wavefunction products can be 
replaced by their ensemble averages, which by Eq. (41) 
makes the internal propagators effectively momentum- 
conserving. Now the analogy to the Gross-Neveu theory 
is complete. 

,/From this analogy we learn that the one-loop RG flow 
computed by two of us 9 is indeed the exact answer in the 
large-g limit, and that the phase transition it indicates 
is real. However, this analogy is even more fruitful in 
providing a way to compute physical quantities in the 
strong-coupling phase 10 , to a description of which we now 
turn. 



B. The Effective Theory of the Strong-Coupling 
Phase 



The strong-coupling phase occurs for u m < u* m for 
any m. We will consider the simplest situation, when 
just one angular momentum channel undergoes an insta- 
bility. This is not an artificial restriction for two reasons: 
Firstly, as seen in the RG treatment, the instabilities in 
different channels are independent. Secondly, it is highly 
unlikely that more than one Landau parameter is close 
to an instability in a generic system. Furthermore, we 
will restrict our attention to charge-channel instabilities 
in the spinful case, in which the spin index is a pas- 
sive spectator. This may turn out to be the experimen- 
tally relevant type of instability; looking at the Landau 
parameters 42 for the clean bulk 2DEG at r s = 5 (Eq. 
(21)) we see that the m = 2 charge channel is the one 
closest to an instability. Instabilities in the spin chan- 
nel are extremely interesting, but in treating them one 
must necessarily also account for the exchange coupling 
J. The resulting theory thus depends on two parameters 
(Z m ,J), and will be investigated in detail in a future 
publication. 

Thus, our model Landau interaction function for this 
section is 



with no sum over m. All the thermodynamic properties 
of the system can be obtained from the partition function 
Z, which can be converted into a path integral in imagi- 
nary time by performing the usual time-slicing according 
to 



Z = Tr(e-< 3H> ) = Tr Yim f[ e -^ H ^ 



(75) 



where (3 = 1/T is the inverse temperature. One then in- 
serts a complete set of fermionic coherent states between 
each factor of e~^ H l N ', and writes the matrix elements 
in terms of the coherent state labels. In the process, 
the fermionic operators c, get replaced by Grassmann 
numbers, which we will represent by ip, The partition 
function can now be written as 



DtpDipe' 



(76) 



where 



0/2 



S = J dt[ + 

-13/2 V " 

_ \u^a ^ cos m{e _ 0/ )nk(t) „ k , {t) \ . (77) 

k,k' / 

Here nu_(t) — tpu.(t)tpk(t) and we have made explicit the 
fact that u m < 0. We factorize the interaction term in 
one time slice (of thickness St = /3/jV): 

exp [St— Xkk' \ u m\ cos(m0 - m0')rik7Jk'] = 



J dcr exp— St A 



..2\0V_ 
1 2\u m \ 



j^iik(ffi cos m6 + a 2 sin m6) 

k 



J dcr exp —g 2 StA 



v 2| Mm | 



where er = (a 1 , er 2 ) has two components, as does M: 
(M 1 , M 2 ) a p = EkC( k )0/3( k ) (cos m6 , sin m6). Note 
that this factorization with a real action for er can be 
carried out only for u m < 0, which is indeed the case 
here. The factors of g are chosen for later convenience. 

The fermionic action is now quadratic, and the 
fermions can be integrated out to obtain an effective ac- 
tion for er 



Seff - -Tr In [(iw - e Q )I - gAcr ■ M] + ,g 2 A J 



dt- 



2|u m | 

(78) 



u m cos m(9 — 9') 



(74) 



where I is the unit matrix and er is now a function of 
imaginary time t. 

To make further progress we expand S e f / in powers of 
er, which can be expressed graphically by the set of "ring 
diagrams" shown in Fig. 3. 

Evaluating the diagram with two external er legs, we 
find the following quadratic contribution from the Tr In 
in the static limit (when fluctuations of er in imaginary 
time are ignored): 
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FIG. 3. The Feynman diagrams in the graphical expansion 
of the Trace Log. Each diagram has one fermion loop with 
different numbers of external cr legs. 



a/3 P k,k' 

((7i cos to(9 + (72 sin m#) (<7i cos m& + sin m#') (79) 



where Np(a) is the Fermi occupation of the single- 
particle state a. 

This quantity can be recognized as the one that ap- 
pears in the internal summation of the leading large- N 
one loop bubble, Eq. (40). As mentioned before (and 
shown in Appendix A) this quantity is self-averaging; 
its disorder-average dominates its fluctuations by g. We 
carry out this self-averaging using RMT to obtain the 
average quadratic part of the action for static cr 



So = g< 



dt- 



M 2 A 



1 



In 2 



(80) 



where the bar represents the fact that we are consider- 
ing only the (self-) disorder-averaged part of So, and the 
In 2 = l/|it£j arises from the summation over energies in 
Eq. (79). For the spinful case this number would be twice 
as large. It is clear that for \u m \ < |tt£j, cr = Q is a stable 
solution, while cr will break symmetry for \u m \ > |w^J. 
In both cases let 



1 



1 



(81) 



denote the distance from criticality. 

Will this symmetry-breaking, present in the "bare" 
theory, survive fluctuations? The answer to this hinges 
on the fact that a factor g 2 appears at all higher loops, 
i.e., in front of the entire self-averaged action for static 
cr. The 4 t,l -order term suffices to clarify the issue. The 
generic 4 th -oidei term is shown in Figure 4, and is repre- 
sented by the following expression (all factors except for 
g have been suppressed): 



4 4 M a pMa~M~sMsa 

9 G ~3 



(82) 



where e is a generic energy denominator. 



FIG. 4. The Feynman diagram corresponding to the 
fourth-order term in the effective action. 



We now replace the product of wavefunction sitting 
inside the M's by their averages, giving the usual argu- 
ments about self-averaging. Just as the average of four 
(/>'s had many terms corresponding to different contrac- 
tions of indices (Eqn. (34)), there are many terms here 
as well. In the leading contribution, the <fip from M a p 
is contracted with the (jf^ from Mg 7 and so on around 
the loop. This forces the same k index to flow around 
the loop, generating the factor X^kVs 4 = Vff 3 - There 
are 4 energy summations (going up to gA) and three fac- 
tors of energy in the denominator (each typically of order 
gA), leading to an extra factor of g. The overall result is 
therefore of order g 2 a 4 . 

It can be shown along very similar lines that the factor 
of g 2 can be extracted from the entire self-averaged effec- 
tive action in the static limit. Thus g 2 plays the role of 
1/Ti (or the inverse temperature), and controls the size of 
fluctuations around the saddle-point. This is what ren- 
ders the large- TV approximation realistic for ballistic dots 
with large dimensionless conductance. 

Armed with the factor of g 2 multiplying the entire 
static self-averaged effective action, one can see that the 
saddle-point approximation becomes exact in the large- 
g limit, and therefore our conclusion that symmetry- 
breaking takes place for \u m \ > |u^J is in fact correct 
in this limit. 

To see what the order parameter corresponds to, re- 
call that in terms of the original fermions ((01,(72)) = 
X)k( n k)(cosm#,sinTO#), where () denotes a quantum- 
mechanical (but not disorder-ensemble) averaging. This 
shows that some values of k are preferentially occupied 
compared to others. Therefore this can be identified 
as the mesoscopic, disordered, version of the Pomer- 
anchuk shape transition of the Fermi surface 27,26 . There 
has been a revival of interest in the bulk transition 
recently 43,44 . Due to the disorder, k is not a good quan- 
tum number, and will suffer quantum fluctuations, 
unlike the clean bulk case. The quantum-mechanical av- 
erage (rik) is similarly not restricted to be or 1. The 
Fermi "surface" will not be sharp, but will be smeared 
out by the disorder. However, the angular average corre- 
sponding to cr still provides an unambiguous measure of 
symmetry-breaking . 
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C. The Connection of Large- TV to the Hartree 
Approximation 

The static large- N saddle-point solution is the same 
as standard Hartree mean field theory (the Fock term 
is down by 1/g because of the form of the Fermi liquid 
interaction). This can be seen very easily by noting that 
the saddle point satisfies the self-consistency condition 

I ( - Tr In [ (iw - e a )l - 9 Aa • M] + g*A J * Jgl) 

(83) 

The above equation in fact minimizes the expression 
for j3 times the free energy, or just the ground state en- 
ergy as T —* 0. Note that it has two parts, one from the 
fcrmionic Tr In and one from the er field energy, ^Tj ^ . 

Using the Feynman-Hellman theorem, which equates 
the er derivative of the expectation value of the fermion 
hamiltonian to the expectation value of the derivative of 
the fermion hamiltonian, we have 

T&^-Tr In[(iw - e a )l - g A<r • M]) = 
gA(SD(cr)\ J2k n ^( cosm0i + smm6j)\SD(a)) (84) 

where we have introduced the Slater Determinantal state 
\SD((t)), which is obtained by diagonalizing the Hamil- 
tonian with a static er, and filling up the lowest g/2 levels. 
This leads to the self-consistency condition 

u m (SD(cr)\ ^ nk(cos77j#i + sin mdj)\SD{a)) = gcr 

k 

(85) 

This is exactly the condition one would get by decou- 
pling the interaction in a Hartree approximation, which 
indicates that at the saddle-point level, the large- TV the- 
ory is the same as the Hartree approximation. However, 
the large- N nature of the theory allows one to justify the 
saddle-point, and to systematically calculate corrections 
to it. 



V. PHASE DIAGRAM AND PHYSICAL 
SIGNATURES 



1/L, and their ratio g ~ L. Thus the limit g — > oo, cor- 
responds L — > oo, i.e., a dot of infinite size. In such a 
system Et — > and everything we discuss (states sensi- 
tive to boundary conditions and a phase transition that 
precedes that of the clean bulk system) occurs within 
an energy window of zero thickness and is thus invisible. 
This is to be expected because in this limit we are really 
describing an infinite ballistic system, t 

We shall refer to the transition that occurs within Et 
in the g — > oo limit as the mesoscopic Pomeranchuk tran- 
sition. With the above words of explanation, the reader 
should not have any trouble reconciling the juxtaposi- 
tion of the words mesoscopic and phase transition in one 
phrase. 

While the g — > oo quantum dot is not directly ac- 
cessible, we are interested in it because several nontriv- 
ial features of the phase transition it undergoes mani- 
fest themselves even for 1/g > 0. This result is based 
on recent investigations into quantum phase transitions, 
which show that while the T = transition (as a func- 
tion of some coupling) occurs at an isolated point that 
is also physically inaccessible, it controls the physical 
behavior in a fan-shaped quantum critical regime 31 in 
the temperature-coupling constant plane. In our prob- 
lem 1/g plays the role of T since g 2 multiplies the ac- 
tion. The situation is depicted in Figure 4. The dot- 
ted lines bounding the quantum critical regime indicate 
the sharp crossover which replaces the transition. There 
are however some differences from the usual quantum 
critical descriptions because g not only appears outside 
the action and plays the role of inverse temperature, 
but also appears within S e ff via subleading, disorder- 
realization-specific terms which are especially important 
in the strong-coupling phase. 

In this section we will explore the three regimes: weak- 
coupling, strong-coupling, and quantum-critical as well 
as crossovers between them. The physics of the weak- 
coupling regime is controlled by the Universal Hamilto- 
nian. However, we will point out that even in a system 
belonging to this regime one can access the physics of 
the quantum-critical regime by making a finite frequency 
measurement. Thus even dots in this region can be used 
to confirm the existence of the phase transition we have 
found at u*. 

We first turn to the strong-coupling and quantum- 
critical regimes which are new to our work. 



We now turn to qualitative and quantitative conse- 
quences of the large-g solution. Our discussion will use 
the language of phase transitions - collective variables be- 
coming gapless and order parameters turning on at crit- 
icality - but there will be inevitable modifications due 
the fact that we are primarily interested in a mesocopic 
system with discrete energy levels. 

Let us begin by reiterating that a genuine phase transi- 
tion occurs only in the limit g — > oo. This is a mathemati- 
cal idealization that allows one to study and classify dots, 
but is inaccessible in practice because A ~ 1/L 2 ,E T ~ 



'This idealization was also implied in RG analysis when it 
gave a zero of the /3 function at u* . In order to get singular 
behavior out of a zero of a /3-function, we need to be able 
to carry out the RG transformation an indefinite number of 
times. This means Et must contain an infinite number of 
states. In a real dot the flow ends when we come to within 
A of the Fermi surface. Thus the critical state with r — 
is characterized by the fact that a has gap of order A as 
compared to the generic value of gA. 
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the quadratic part of the effective action for slow fluctu- 
ations of cr. We find 



critical fan 



broken 



u* symmetric 



1/fl 



FIG. 5. The generic phase diagram for a second-order 
quantum phase transition. The horizontal axis represents the 
coupling constant which can be tuned to take one across the 
transition. The vertical axis is usually the temperature in 
bulk quantum systems, but is 1/g here, since in our system 
one of the roles played by g is that of the inverse temperature. 



A. The Strong-Coupling Regime 

As shown by Eq. (80) cr acquires a nonzero expecta- 
tion value if u m < u* m or r < 0. Since cr can be thought 
of as a two-dimensional vector order parameter, and the 
self- averaging part of S e ff has rotational symmetry, one 
expects a Mexican Hat (a circle of exactly degenerate 
minima |<x| = <tq). For g = oo the system can still break 
symmetry and pick a particular value of cr, since there are 
no kinetic term for cr in that limit (see Eq. (87) below). 
However, for 1/j ^ one may expect that cr will de- 
localize immediately around the Mexican Hat, restoring 
symmetry, since 1 / g 2 plays the role of T or h (depend- 
ing on how one wants to interpret the path integral). In 
other words one expects the zero-dimensional analog of 
the Mermin- Wagner theorem 49 . In terms of the parame- 
terization of cr in cylindrical coordinates 



& = fa X) 



(86) 



we expect the cr variable to uniformly populate all values 
of the angle x- 

This does not actually happen, thanks to disorder- 
realization-specific terms in S e tf, which are subleading 
in 1 1 g to the terms considered so far and explicitly break 
rotational symmetry in cr space. There are many sources 
for such 1/g terms. The coefficient of the quadratic term 
coming from the Tr In, which has a leading average part 
which produces the In 2, also has sample-to-samplc fluc- 
tuations which are down by 1/g (as shown in Appendix 
A). In fact, a term first- order in cr is also present for a 
system undergoing an instability in an even angular mo- 
mentum channel m. These subleading terms render the 
valley of the Mexican Hat nondegenerate in the case of 
even m, and leave behind a two-fold degeneracy in the 
case of to odd. 

Will these sample specific terms localize cr in the angu- 
lar direction? To resolve this issue we need to understand 
the low-energy dynamics of cr. To this end we compute 



So = g 2 



dt- 



£W)i a /(«) 

(87) 



where the bar represents the fact that we are consider- 
ing only the (self-) disorder-averaged part of Sq. The 
frequency dependent function f(tu) is 



2wtan~ 1 f 



Aln 



Let us note that 



(4^ J +A 2 ) 



(w 2 +A 2 ) 



Aln 



^tan" 1 ^ 



1 4- — 

' A 2 



f(oj) = uj 2 /AA 
= ttM/4 



oj < A 
A < oj < E T 



(89) 
(90) 



Thus we expect Hamiltonian dynamics at very low ener- 
gies, with a non-universal sample specific dependance on 
the lower cutoff A. At high energies we have overdamped 
dynamics, representing the effect of Landau damping on 
the cr variable. This latter behavior is universal and self- 
averaging. * 

At very low energies described by hamiltonian dynam- 
ics Eqs. (89,) the momentum conjugate to cr is 



Per = g&/2A. 



(91) 



Only the angular dynamics is important in the 
symmetry-broken phase since the radial fluctuations of 
cr are suppressed by the Mexican Hat potential of order 
g 2 A. The Hamiltonian for small oscillations around the 
global minimum of the explicit symmetry-broken effec- 
tive action looks like 



H osc = A h gAV(\) 

9 



(92) 



where we have denoted the angular part of the canonical 
momentum as L. The spread of the ground state wave 
function (of this oscillator) can be estimated as 5\ — 
1/y/g- Thus at large g the support of cr is localized to 
the minimum of the potential energy term. Our analysis 
is internally consistent: The form of f(u>) we used to 



'The reader might ask what permits us to keep some terms 
which are 1/g down compared to the self-averaged static part 
of S e ff, and discard others. The answer is that we keep only 
the lowest order terms of every type. Thus, the realization- 
specific terms which break the rotational symmetry of the 
self-averaged part of S e ff are kept, while 1/g corrections to 
the self-averaged part are discarded. Similarly, though the 
kinetic term of Eq. (87) is order 1/g down compared to the 
self-averaged part, it must be kept since this is the lowest 
order dynamics. 
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obtain oscillator description (Eqn. (89)) is valid out to 
u ~ A, the oscillator ground state energy. § 

Thus symmetry-breaking will occur in the strong- 
coupling regime for large enough g when disorder- 
realization-specific terms are included. The result is a 
consequence of spontaneous symmetry-breaking (requir- 
ing strong interactions) and explicit symmetry-breaking 
(from disorder). 

These considerations apply uniformly for m even and 
odd. However, there are significant differences in other 
features and we will treat them separately. 

1. The Case of m Even 

Let us see how to find the effective potential with the 
realization-specific corrections more explicitly. We need 
to find the "potential landscape" of er, which we define 
for t-independent er as 

er 2 1 

Veff(v) = g 2 A- - -Trln[(iw - e a )S afj - .gAer ■ M afj ] 

£\ u m\ P 

(93) 

where (3 = 1/T is the inverse temperature not to be con- 
fused with the label for noninteracting disorder eigen- 
states. We thus need to map out the Tr In for time- 
independent er. In the limit of zero temperature, the 
quantity Tr ln(iu> — H Fl cr)/P is just — £p(er), where the 
fermionic Hamiltonian is 

Hf.ct = e " c l c " - 9 Act ■ M a/3clcf3 (94) 

a afi 

and £f{&) is its ground state energy, obtained by diag- 
onalizing the hamiltonian Hf.ct and filling up the lowest 
states with the requisite number of fcrmions. Let us note 
that for large g, since fluctuations in er are small, we 
can approximate the total energy of the interacting sys- 
tem with N particles £n as the minimum of the effective 
potential of Eq. (93), 

( 2 A _2 \ 
£ f {<t) + V— (95) 

where Miner indicates that the argument should be min- 
imized over all er. Note that the value of er at the min- 
imum depends on the realization as well as the number 
of particles N. 

For small er we can find £p by standard perturbation 
theory. To first order we find the following energy cor- 
rection 

£ { f\v,X)= -gAtr-^M aa (96) 

aocc 

= -gAo £ 0* (k)0 a (k) cos(m0 - X ) (97) 

k,aocc 



''The low frequency approximation f(oS) ~ oj 2 /(4A) tracks 
the exact expression to with ten percent in this range. 



where the sum is over occupied states only, and a, x are 
the magnitude and angle of er respectively. From the fact 
that </>* (k) = 4>a{— k) it follows that 

M aP = (-l) m M; a (98) 

a specific instance of which is 

M aa = (-l) m M QQ (99) 

This in turn shows that M QQ must vanish for odd m. 
However, for even m there is no reason for this to vanish 
for a particular disorder realization (though its disorder- 
average does vanish). One can also see that this term 
produces a dependence on the angle x of er, which breaks 
the degeneracy of states in the Mexican Hat. To get an 
idea of the magnitude of this realization-dependent term, 
we square it and average over disorder. It is easy to verify 
that the disorder ensemble average (denoted by <C^>)is 
given by 

< (4 X) ) 2 >- 9 2 A 2 (100) 

meaning that the typical size of this correction is of order 
g, one power of g down from the self-averaging part of £ . 

One can also analytically get the statistical correla- 
tions in the shape of the first-order contribution to the 
potential by computing 

« 4 1) (Xi)4 1) (X2) »^ .9 2 A 2 cos( X i - X2) (101) 

The picture we get for even m is then the following: 
Each order in perturbation theory contributes a fairly 
smooth (on average) realization-dependent term which 
depends on x, thus breaking the symmetry of the g = oo 
Mexican Hat. If the system is not too deep in the strong- 
coupling phase we can expect that only a few orders 
of perturbation theory contribute, and that the angle- 
dependent potential in the Mexican Hat is fairly smooth, 
and gcncrically has a single minimum. The most impor- 
tant effect of the realization-dependent terms is to allow 
for symmetry-breaking even for finite g. Therefore, both 
disorder and interactions are crucial to this symmetry- 
breaking: Interactions produce the Mexican Hat potential 
causing spontaneous symmetry-breaking at g = oo, while 
explicit sample specific symmetry breaking terms prevent 
the restoration of symmetry by quantum fluctuations even 
at finite g. 

At T = one can simply evaluate the fermionic ground 
state energy £p{cr) corresponding to each static er nu- 
merically. We show the resulting effective potential of 
Eq. (93) in Fig. 6 for the illustrative case m = 2. As can 
be seen, there is an approximate Mexican Hat circle of 
degenerate minima. In Fig. 7 we zoom in on this circle 
to see that the degeneracy is broken by sample-specific 
disorder and that the fluctuations are indeed of order 1/g 
relative to the self-averaging term. Clearly, the numer- 
ical work corroborates the analytical picture developed 
above. 

One can now ask how various measurable quantities 
behave in the strong-coupling phase. There are three 
such quantities of interest. 



16 



c 0.8 
o 




6 




FIG. 6. The effective potential as a function of a for m = 2, 
with g = 20 and u m = — 1.2 (left panel) and u m — — 1.7 (right 
panel). The symmetry-breaking in the second case is clear. 
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FIG. 7. The effective potential in the Mexican Hat in units 
of A plotted as a function of the angle \ f° r m — 2, at g = 20 
and Um =1.7. Note the single nondegenerate minimum and 
the scale of sample specific disorder, which is of order 1/g 
relative to the average. 



FIG. 8. The peak-spacing distribution for m = 2, at g = 20 
and w m = 1.7. Note the incidence of negative values of the 
peak-spacing, which is impossible in a noninteracting model. 

• The addition spectrum, which measures the differ- 
ence in ground state energies Ai(N) = £n+i ~ £n 
where the total energy of the system with N parti- 
cles £m is the sum of the fcrmionic energy and the 
energy of the bosonic field er (Eq. (95)). We have 
not included the charging energy UqN 2 /2. The 
Coulomb Blockade peak spacing is then given by 
A 2 (N) = Ai(JV + l)- Ai(iV). 

• The peak-height distribution. 

• The persistent current in response to an external 
flux. 

Since at any particular er the electronic excitations 
are controlled by a single-particle Hamiltonian, and er is 
well-localized, one might naively imagine that the addi- 
tion spectrum is the same as the level-spacing spectrum, 
which is the Wigner-Dyson distribution with an enhanced 
A. However, this ignores possible shifts of er upon adding 
a particle. Most of the time the value of er shifts only a lit- 
tle upon adding a particle. However, every now and then 
a minimum which was far away in er-space becomes lower 
upon adding a particle. One can compute the peak spac- 
ing from A 2 (N) = Ai (iV + 1) - Ai (iV). In Figure 8 we 
show the distribution of A2 calculated numerically by the 
method explained above Eq. (94). As can be seen, nega- 
tive values of A 2 are possible because of the rare events in 
which er changes discontinuously upon adding a particle, 
which makes the distribution broader and more symmet- 
ric than the level-spacing distribution. (Note once again 
that the negative answer results only upon subtracting 
out the charging energy Uq.) In contrast negative values 
of A2 are impossible in a noninteracting model. 

Now consider the peak heights. In the noninteract- 
ing case, the distribution has been computed in RMT 
by Jalabert, Stone, and Alhassid 39 , and measured ex- 
perimentally in the weak-coupling regime at very low 
temperatures 40 . There are characteristic differences 39 
between the peak-height distributions for the GOE 
(time-reversal unbroken) and GUE( time-reversal bro- 
ken) which have been measured 40 , and will play a role 
in the following subsection. There are discrepancies be- 
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FIG. 9. The ensemble-averaged ground state energy as a 
function of external flux for m = 1 (solid line) and m = 2 
(dashed line), at g = 20 and u m = 1.7. Note the fact that the 
ground state energy increases for m = 2, which thus shows 
diamagnetic behavior, while m = 1 shows paramagnetic be- 
havior. 

tween theory and experiment at higher temperatures 50 
which we will address qualitatively in a later subsection. 
As long as er changes by small amounts when particles are 
added the strongly-coupled system behaves much like the 
noninteracting system. However, when a large change of 
er is involved, all the electrons must change their state 
to adjust to the new a. This means the transition ma- 
trix element the fermion creation operator between the 
two ground states will be small, and we expect the peak 
height to be considerably smaller.** This effect produces 
a correlation between smaller than average peak spacing 
and small peak height. 

Finally, let us briefly consider the response of the sys- 
tem to an external flux, leaving the details to Part II 11 
(a short report of this work can be found in ref. 12 ). This 
response is none other than the persistent current 

dF 

Ipers — q , (102) 

where F is the Free energy (the ground state energy at 
T = 0). This response is primarily orbital in GaAs due 
to the tiny coupling of the external flux to the spins. 

Figure 9 shows the ground state energy as a function 
of external flux averaged over many disorder realizations. 
The noteworthy point is that the average persistent cur- 
rent is diamagnetic. To see why this is special, a brief 
digression into what is known about persistent currents 
is necessary. In a mesoscopic ring penetrated by a flux, 
the ground state energy has to be periodic in the flux, 
since an integer number of flux quanta can be gauged 
away. So 



** Changing er is the same as changing a parameter in a non- 
interacting random Hamiltonian, a theme which has been ex- 
plored in the literature 61 . Borrowing the result for the overlap 
of determinants 62 , for ground stated differing by Scr we expect 
the peak height to scale as e -9 ln 9(«°") 2 . 




= 2i sin(27T0/0 o ) + h sin(47rt/>/(fo) H 



(103) 

where 4>q = h/e is the flux quantum. 

The noninteracting problem is relatively well- 
understood 2 ' 51 . Only the even moments I 2n survive 
disorder-averaging (this result holds for the interacting 
case as well). The typical, fluctuating values of the 
Fourier coefficients are (for small n) I n ,typ ~ Et/4>o while 
the average is (hn) ~ A/(/> . 

Interactions, when included in renormalized first-order 
perturbation theory 32 , produce 

(h) ~ L i *~T~ ( 104 ) 

where fj,* (typically < 1) is the dimensionless Cooper- 
channel interaction at low energies. The conclusion is 
that interactions enhance the average persistent current, 
but that if /i* > (I2) should be paramagnetic, while if 
fx* < it should be diamagnetic. Much numerical work 
has followed since 52 , mostly on one-dimensional rings, 
confirming that interactions do indeed enhance (I2) in 
the spinful case. 

Experiments 8 reveal striking discrepancies with the 
above predictions. The predicted value 32 , while of the 
right order-of-magnitude, is still smaller than experi- 
ment. More disturbingly, the sign is inconsistent. Even 
materials that show no sign of superconductivity (im- 
plying that n* > 0) show 8 a diamagnetic (12). Many 
other explanations (e.g. ref. 53 ' 54 ) have been proposed to 
account for these observations, but questions about the 
sign and the magnitude of persistent currents in non- 
superconducting materials remain open 33 . 

In this context it is striking to obtain, as we do for m 
even, a diamagnetic persistent current in a model without 
superconductivity. 

All the signatures described above (expect for the dia- 
magnetic persistent current) also hold for the case of m 
odd. However, there are further interesting effects which 
we describe below. 



2. The Case of m Odd 

Let us now turn to odd m, where the first-order con- 
tribution in er to the effective potential is missing. In 
fact, one can show a much more general result, that the 
ground state energy to all orders £{x) &t a particular 
angle \ with a = \tr\ fixed satisfies 

£(x) = £(x + tt) (105) 

To show this result we write the second-quantized 
fcrmionic Hamiltonian for the illustrative case of to = 1 
as 

Hf(x) = ^£ a cic a ~gAa ^ 4c/3C( k )0/3(k) cos(6> - x) 

(106) 
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and its first-quantized version: 



ha/j = 8 aP e a - gAcr^0*(k)^(k)cos(6l - x). (107) 

k 

Using Eqs. (32,98) along with 0(-k) = 0(k) ± tt one 
can easily show that 

h*(x)=h(x + n) (108) 

Since the Hamiltonian is hermitian at every x its eigen- 
values are real, and thus 

Hx)^n(x) = e«(x)*n(x) 
h*( X )K(x) = e n (x)K(x) 
^h{x + K)K{x)=e n { X )K{x) (109) 

which shows that all the single-particle energies of h(x) 
are shared by h(x + tt), so naturally the ground state en- 
ergies are the same as well. This implies that the energy 
as a function of angle x generically has two exactly de- 
generate minima. The complex conjugate relationship 
between the single-particle wavefunctions shows that the 
two minima are connected by a time-reversal transfor- 
mation T. This is easy to understand for m = 1. In the 
bulk an m = 1 distortion is just a shift of the Fermi sur- 
face. Since the underlying Hamiltonian is time-reversal 
symmetric, any distortion should have the same energy 
as the time-reversed distortion. 

It is also important to understand that an m = 1 dis- 
tortion in the quantum dot does not represent a current 
flowing in a particular direction, as it would in the bulk: 
This is impossible in the steady state in a closed dot. 
Rather, the m = 1 distortion describes a persistent cir- 
culating current within the dot. One can understand this 
by referring to Eq. (29), which shows that the state we 
label by k in the dot has a vanishing wavefunction at 
the boundary (assuming a hard wall boundary condition 
for simplicity), and therefore the current carried by it 
vanishes at the boundary of the dot. 

Once again, one can numerically corroborate all these 
deductions, as shown in Figs. 10 and 11, which show the 
approximate Mexican Hat and the double-degeneracy in 
it for the illustrative case of m = 1. 

Similarly, the peak-spacing distribution can also be 
numerically obtained, as shown in Fig. 12. 

Thus, the 0(2) symmetry of the Mexican Hat in the 
g = oo limit is broken down to a two-fold degeneracy 
for any finite g. This leads to a very interesting physics 
which depends sensitively on the coupling to the leads. 

Isolated Dot 

The isolated dot has two exactly degenerate minima in 
the effective potential. However, the low-energy dynam- 
ics of <r will induce a "tunneling" term between the two 
minima. Recall that the "kinetic" part of the effective 
action at very low energies was approximately 

S e ff,kin = 9 J dr (^^j (HO) 




FIG. 10. The effective potential as a function of cr for 
m = 1, with g = 20 and u m = —1.2 (left panel) and 
Um = —1.7 (right panel). 
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FIG. 11. The effective potential in the Mexican Hat in 
units of A plotted as a function of the angle x f° r the case 
m = 1, at g = 20 and u m = 1.7. Note the exact two- fold 
degeneracy of the effective potential minima, and the scale of 
sample specific disorder, which is of order 1/g relative to the 
average. 
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FIG. 12. The peak-spacing distribution for m = 1, at 
g = 20 and u m =1.7. Note the incidence of negative values 
of the peak-spacing, which is impossible in a noninteracting 
model. 
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The tunneling amplitude typically goes as 

A = e ~I dx ^ 2mV{x) (111) 

where the integral is over the classically forbidden region. 
In this problem the mass, read off from the kinetic term, 
is of order g, and so is the potential in the angular di- 
rection. Thus the tunneling amplitude must be of order 

Thus the true eigenstates of the system will be the 
symmetric superposition, which is the ground state, and 
the antisymmetric superposition which is the excited 
state. The splitting between them will be of order 
e~fA.The ground state at T = will not have any 
symmetry-breaking . 

However this state of affairs describes a very narrow 
low temperature range: For A>T> e~ 9 A the system 
will be in an (essentially) equal incoherent mixture of the 
symmetric and antisymmetric states, or equivalently, by 
change of basis, in an equal incoherent mixture of the two 
minima of the effective potential. Each of these minima 
has T-broken dynamics and GUE statistics 39 ' 40 even in 
the absence of an external magnetic field. Since the two 
minima have complex conjugate wavefunctions, the inco- 
herent average will wipe out the the persistent currents 
which change sign under complex conjugation, but not 
the GUE peak-height statistics 39 which are immune to 
complex conjugation. 

Another way, besides increasing T to experimentally 
observe this degeneracy is to explicitly break time- 
reversal symmetry by introducing a small external mag- 
netic flux <j>. The response of the system to such a flux 
will be linear instead of quadratic, as in the time-reversal 
invariant cases, here are some specifics. 

The flux is introduced within the theory by replacing 
the non-interacting part of the Hamiltonian by one drawn 
from an ensemble of crossover Hamiltonians 14 . Including 
the external flux, the first quantized Hamiltonian at a 
particular angle \ is (from Eq. (106)) 

h = h(x) + ^h A (112) 
<Po 

where C is a constant of order unity that depends 
on the shape of the dot and the nature of boundary 
scattering 21 , and Iia is an appropriately normalized, ran- 
dom, pure imaginary, antisymmetric matrix in any or- 
thogonal basis 14 with typical matrix elements of size 
A^/g. It is easy to show that the two minima of V e ff 
move in opposite directions to first order in <f>. It is con- 
venient to use the exact noninteracting eigenstates as a 
preferred orthogonal basis. Using the notation of Eq. 
(109), the change in energy to first order is 

S£ F (X) = J2 £(*n,a(x))*>lA,a|8* n ,/j(x) (113) 

n occ a/3 

where $ n ,a(x)) is the projection of the n th exact eigen- 
function of Eq. (112) on the noninteracting disorder 
eigenvector a. Noting from Eq. (109) that (3> n ,a(x))* = 



^n,a(x + 7r )j an d using the antisymmetry of Ha in the 
basis of exact noninteracting eigenstates, it follows that 

S£ F (x) = -S£ F ( X + 7r) (114) 

One can estimate the magnitude of 8£p by squaring 
and averaging over the ensemble to be 8£p ~ (gA)<p/<f>o- 
When this energy exceeds the splitting between the sym- 
metric and antisymmetric states, the system will lo- 
calize cr in the lower minimum of V e ff, thus breaking 
time-reversal symmetry completely. The flux needed to 
achieve this at T = is known as the crossover flux 

Across ^ <Po^ 9 /g (115) 

which is exponentially small in g. The effective schematic 
hamiltonian in the subspace of the two minima, labeled 
+ and — is of the form 

< 116 > 

where b denotes the external magnetic field and Ae~ 9 
the tunneling amplitude between the two minima related 
by time-reversal. 

This response to an applied flux should be contrasted 
with the noninteracting case, when the crossover flux (for 
single-particle properties) is of the order of </>o / ^fg. The 
persistent current in the m odd case is paramagnetic 12 , 
as can be seen from Fig. (9). The ground state energy 
decreases linearly for small 0, a consequence of the first 
order contribution. 

The response of the time-reversal broken system man- 
ifests itself in the shift in the zero-bias conductance peak 
position as a function of external flux. Since the peak 
position reflects the difference Ai = £n+i — £n, and all 
the single-particle energies change linearly with flux (for 
fio/y/g S> ^ Across), the shift in the peak position 
must also be linear in the external flux. This contrasts 
with the noninteracting case where the shift is quadratic 
in the external flux for (j>o/y/g ^> <t>- The numerical re- 
sults are shown in Fig 13. 

For A 3> T 3> e~ 9 A one can carry out a linear- 
response magnetization measurement, about which more 
will be said in Part II 11 . The two degenerate states will 
lead to a 1/T Curie- like linear response. Let us now turn 
to what happens when the quantum dot is weakly cou- 
pled to the leads. 

Dot Coupled to Leads: Mapping to 
the Caldcira-Leggctt and Kondo Problems for m odd 

We have already seen that at very low energies the 
dynamics of cr in an isolated dot becomes Hamiltonian, 
due to the fact that there are no particle-hole excitations 
below a lowest energy. This changes if the dot is coupled 
to leads, since the dot levels will hybridize with the con- 
tinuum of levels in the leads. This leads to dissipative 
dynamics for cr at arbitrarily low energies. Since in the 
case of m odd we know that there are two degenerate 
minima of the effective potential, we are naturally led to 
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FIG. 13. The peak position Ai = £n+i — £n as a function 
of external flux in the range < C<f> < <j>o/^/g for m — 1, at 
g = 50 and u m = 1.7. Note the fact that the peak position 
is linear in the flux for small values of flux in contrast to the 
noninteracting result (dashed line) which is quadratic. 



the Caldeira-Leggett problem 34 , which concerns a parti- 
cle in a degenerate double- well potential interacting with 
a "bath" of harmonic oscillators which model the dissipa- 
tion. Applying the solution 35,36 of the Caldeira-Leggett 
model to our problem, we can predict that there will be 
two phases depending on the strength of the coupling to 
the leads: One with spontaneously broken time-reversal 
symmetry where the cr field is trapped in one minimum, 
and one where it can successfully tunnel back and forth 
between the two minima and time-reversal symmetry is 
unbroken. 

Let us recall some of the principal results of Caldeira 
and Leggett 34 . Their motivation was to study systems 
with a macroscopic quantum mechanical degree of free- 
dom (such as the flux through a superconducting loop or 
cr in our case) coupled to a dissipative environment. If 
the particle is initially set up to be in the lowest state 
of a single well, in the absence of dissipation the parti- 
cle will eventually tunnel between the minima, with the 
ground state being a symmetric superposition of the low- 
est states of the two wells. In the presence of dissipation 
tunnelling becomes sensitive to the spectral density of the 
harmonic oscillator bath at low frequencies. For the case 
of ohmic dissipation, Caldeira and Leggett concluded 34 
that dissipation decreases the tunnelling rate. 

Soon after this, Chakravarty 35 and Bray and Moore 36 
took the results much further. They pointed out a deep 
connection between this problem and the anisotropic 
Kondo problem 55 of electrons interacting with a local- 
ized spin. They reached the important conclusion that 
there are two phases in the Caldeira-Leggett problem. At 
sufficiently weak dissipation the particle is delocalized be- 
tween the two degenerate wells, and this corresponds to 
the Kondo singlet phase of the antiferromagnetic Kondo 
model (Jj_ > 0). At strong dissipation the particle is 
localized in one well, corresponding to the ferromagnetic 
Kondo model (J± < 0), for which the spin-flip terms are 
irrelevant at low energies. 

The relation of this work to ours is made more 
precise if we consider the effective action 56 used by 



Chakravarty 35 for a resistively shunted superconducting 
quantum interference device (SQUID): 



0h 

S eff = S / dt 
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(117) 



where 8 denotes the flux (in units of <Pq/(2tt)), C is the 
capacitance of the SQUID, u{9) is the double-well poten- 
tial with degenerate minima at ±#o, and the dimension- 
less number r\ — h~/2ire 2 R parameterizes the dissipation 
(R is the shunt resistance). The phase transition occurs 35 
(in the limit when the tunnelling amplitude is very small) 
for 4rysin 2 (^o/2) = 1. The effective action for the Kondo 
problem 55 can also be cast in the above formJT 

In our problem, dissipation at very low frequencies 
will enter the effective action through the susceptibil- 
ity, which is the quadratic term in cr in the fermionic 
TrLn. We have seen before (Eqn. (88)), in the case of the 
isolated dot, that the dynamics of cr become dissipative 
for ui » A when it can decay into numerous fermionic 
many-body states. When connected to the leads, this 
Landau damping behaviour persists down to the Fermi 
energy because the dot states are now embedded in the 
continuum and A% h , the spectral densities of the state 
labelled by a above (p for particle) and below (h for hole) 
the chemical potential in the exact basis, evolve from the 
delta function to a Lorentzian of width Y. 



A a {e)=Q{e)Al{e) + Q{-e)A h a {e) 
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r 2 + (e - e a ) 2 

(118) 



Since we are focusing on very low frequencies, the spec- 
tral densities are of order T/A 2 . Remembering that 
| M a p | 2 is typically of order 1/g, and putting in the over- 
all factor of g 2 in the quadratic term, one can easily check 
that the imaginary part of the susceptibility which repre- 
sents dissipation will therefore be of order gT 2 uj. Going 
over to the imaginary time representation, one can show 
that up to an unimportant constant we have the following 
term in the effective action 

, r 2 r dui 
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(119) 



^The dissipation induced by the bath manifests itself as a 
long-range 1/|£ — t'\ 2 interaction in Euclidean time. Exactly 
the same effective action is obtained by integrating out the 
dissipative electrons in the Kondo problem, as was done by 
Anderson, Yuval, and Hamann j5 . At a physical level, the 
particle spends most of the (Euclidean) time in one of the 
wells, occasionally tunnelling between them, an event called 
an instanton. These instantons correspond to spin-flip events 
in the Kondo problem 55 . 
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In addition we also have the yyfcj term from the higher 
energy states. 

Comparing our full effective action (Eq. (87) + Eq. 
(119)) to that of Chakravarty 35 (Eq. (117)), we see that 
the angle x of <x corresponds to the "particle" degree of 
freedom 0, and that gT 2 / A 2 w 1 corresponds to large 
dissipation (in the strong coupling phase (a) is of order 
unity). Therefore, despite the fact that the two min- 
ima are exactly degenerate, cr will be localized in a sin- 
gle minimum at large g for T > A/^fg. Even when 
this condition is satisfied the level width is still much 
smaller than the level spacing for large g. Under these 
conditions time-reversal symmetry will be spontaneously 
broken for m odd and zero temperature. As mentioned 
before, time-reversal symmetry breaking can be inferred 
from the statistics of peak heights 39 ' 40 . This symmetry- 
breaking also induces a spontaneous persistent current at 
zero flux, and will be described in Part II 11 . Note how- 
ever that in order to see this effect we need T < Ae~ 9 
since above this T, thermal effects will produce the same 
effect. 

The T-breaking transition is structurally identical to 
the ferromagnetic-antiferromagnetic Kondo transition, 
but is different in character from other Kondo-like states 
proposed in quantum dots 57 ' 58 . Any Kondo-like state has 
to start from a finitely-degenerate set of states forming 
a pseudospin, and "conduction" electrons with a con- 
tinuum of energies coupling to this pseudospin. The two 
main previous proposals are: (i) The degeneracy is in the 
number of particles at the zero-bias conductance peak, 
when the dot has the same free energy to have N or N + 1 
particles, which are the two states of the pseudospin- 1 ;. 
Electrons hopping from the leads "flip" the pseudospin, 
and at low enough T form a Kondo resonance at the 
Fermi energy 57 , (ii) A dot with an odd number of elec- 
trons has a singly-occupied state as the highest occupied 
level with spin-i, which then acts as a standard Kondo 
impurity 58 . In both these examples, one cannot change 
the sign of the Kondo coupling to go from the antifer- 
romagnetic to the ferromagnetic Kondo model. Thus, 
while there is an observable resonance at the Fermi en- 
ergy below the Kondo temperature Tk, there is no phase 
transition (see, however, Ref. 59 ). The new features of our 
state are that the pseudospin is a truly collective variable, 
which represents the Fermi surface distortion of all the 
electrons, and that one can drive the transition between 
the analogs of the antifcrromagnetic and ferromagnetic 
Kondo problem merely by tuning the coupling to the 
leads. In contrast to the g = oo mcsoscopic Pomeranchuk 
transition which broadened into a crossover for any finite 
<?, this is a true quantum phase transition. The finite 
system is able to do this by coupling coherently to the 
infinite reservoir via the leads. 



B. The Quantum Critical Regime 

Now we turn to a description of the quantum criti- 
cal regime, which can be reached from the strong- or the 



weak-coupling regimes by increasing either the probe fre- 
quency,gate voltage, temperature T or l/g. We will focus 
below on increasing l/g. 

The most important feature of this regime is that 
a has a gap of order A whereas in the weak-coupling 
regimes it has a gap of order gA. This in turn leads to 
a broad spectral function for low-energy excited states. 
This could be seen in tunnelling at finite bias 60 (nonlinear 
conductance) . 

1. Connection to Fock-Space Delocalization 

The physics of the quantum critical regime can 
be connected to ideas of Fock-Space localization- 
delocalization 37,38 , which examine the spectral function 
of an excited quasiparticle state in a mesoscopic system 
as interactions are turned on. Assuming no correlation 
between the matrix elements of the interaction, V a ^g, 
these studies concluded that for excitation energy 

e > e* w Axjg/\ng (120) 

the spectral function of the quasiparticle is broadened 
into a Lorentzian, characteristic of the infinite system. 
Below e* the spectral function splits up into a number 
of sharp peaks. We briefly recount the simple arguments 
behind this result (modulo the hit/) 37 as it relates also to 
our work. Consider a single particle of energy e decaying 
into a particle of lower energy and a particle-hole pair by 
Fermi's golden rule. By energy conservation we need only 
consider the energy of the particle-hole pair as indepen- 
dent. The matrix element for the process is V Q( g 7< 5 ~ A/g, 
while the density of states is just the density of particle- 
hole states (or two-body states), which is given by the 
number of ways of partitioning the total energy e among 
three particles ^(e) = £ 2 /2A 3 . Note that this is much 
smaller than the density of many-body states. The total 
width induced by this process is roughly V ~ e 2 /g 2 A. If 
this is much smaller than A one cannot make the con- 
tinuum approximation, and the quasiparticle weight will 
split into a number of sharp peaks. The correct criterion 
is to compare the matrix element between two different 
three-body states \u m \ A/g to the typical spacing between 
such states I/1S3 (e). If the two are comparable, the quasi- 
particle weight gets distributed among many such states. 
This leads to the criterion e* ~ A^fg. To obtain the 
extra lng in Eq. (120) requires a more sophisticated 
analysis 37 , in which Fock-Space is thought of as a lat- 
tice, each site of which corresponds to a single Slater 
determinant. Hopping on the lattice is accomplished by 
the interaction matrix elements V a pfS- If the system is 
localized on this lattice, then the quasiparticle state is 
comprised of at most a few determinants, and the spec- 
tral function will therefore be a sum of sharp peaks. If 
the system is delocalized on this lattice, the quasiparticle 
state has an infinite number of Slater determinants, and 
the spectral function is broadened into a Lorenztian. To 
find the threshold energy for delocalization e* one carries 
out the analog of the Anderson locater expansion. 
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FIG. 14. The lowest-order sell-energy diagram lor the 
lermion due to its coupling to the collective field cr. Straight 
lines represent lermions while wavy dotted lines represent cr. 



Let us see how this gets modified in our approach 
wherein we keep correlations between different V a /3js (all 
of which really arise from a single Landau parameter). 
This gives rise to loop corrections which were neglected 
earlier 37 ' 38 , which lead at the critical point to V a ^& 
which increase as the energy scale is lowered, as discussed 
in Section III. In fact, the Fermi-surface distortion be- 
comes the bosonic dynamical variable cr, which interacts 
with the fermion to give it a decay width. 

To make this more quantitative, let us consider the 
retarded, real-time, inverse propagator of cr (after self- 
averaging the wavefunction sums) in the weak-coupling 
phase 



D ret(u) = i r 7T 
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(121) 



where r\ is a small positive infinitesimal. 

The fermions couple to cr through the term — gAcr ■ 
M (Eq. (78)). The fermion propagator will get dressed 
by this coupling, and acquire a self-energy S ret)Q , the 
simplest diagram for which is shown in Fig. (14). The 
decay width produced by this diagram can be written as 

T a {cj) = - JmE ret , a (w) = g 2 A 2 £ \M a[j \ 2 x 

P 

J &-Im D ret (uj - oj') x ImGp(u>') (122) 



where Im stands for the imaginary part of the function 
that follows. 

In computing the width by the above one loop dia- 
gram, we seem to be following standard diagrammatic 
perturbation theory. There is, however, a subtlety that 
needs to be borne in mind. The er propagator we use is 
the not the bare one based on the last term in the action 
of Eqn. (78), (where it had no oj dependence) but an 
effective one Eqn. (88), obtained upon integrating out the 
fermions. One may ask what the fermion is doing in the 
one-loop graph above if it has already been integrated 
out. Is there some double-counting? In Appendix B we 
will show that there is no double-counting and that the 



above procedure for finding T is legitimate. Appendix B 
also explains how to incorporate symmetry-breaking (a 
nonzero average for er) into the calculation. 

Resuming our calculation of T, we will assume (and 
justify later) that the fermion spectral function is still 
sharp on the scale of a few A so that we can make the 
replacement 



Im(G ret> p(oj)) ~ -ttS(uj - ep) 



(123) 



inside the integral of Eq. (122). We now consider the 
spectral function of the er-propagator. In Eq. (121), for 
A <C oj <C Et, the w-independent part of the sum over 
a/3 produces g 2 /\u* n \, while the dominant w-dependence 
is imaginary. Reading off the imaginary part from Eq. 
(121), we get 



-ioj + C 



where 



and 



C = 2rgA/ir 



(124) 



(125) 



(126) 



is the distance from the critical point. Note that this in- 
volves replacing the discrete sum over particle-hole states 
by a continuum, and will be justified self-consistently in 
the following. The imaginary part of this retarded prop- 
agator is 



Im(D ret (uj)) 



gir to 2 + C 



(127) 



This shows that £ acts like a gap or a mass for cr, and 
when C — > 0, cr should become gapless. Substituting 
Eqs.(123,127) into Eq. (122) we find the decay width to 
be 



r(w) ~ -ln(l +uj 2 /C). 

7T 



(128) 



In the weak-coupling phase at low energies we have 
oj C- Expanding the In and noting that £ ~ g/\u m \, 
we recover the weak-coupling result. At the critical point 
one must replace £ by A. This is the infrared cutoff of 
the theory, signifying the energy below which D ret has 
no imaginary part. This implies that in the quantum 
critical regime, putting the particle on shell (to — e), we 
have 



2A 

r(e) ~ — lne/A. 

7T 



(129) 



Now consider for r(e) in the strong coupling phase. 
We must expand the theory to quadratic order in cr at the 
saddle point. The fermionic hamiltonian at the saddle 
point will take the form 



(130) 
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where we have chosen the direction of explicit (sample 
specific) symmetry breaking to be along the x-axis. The 
fluctuation about this point have a huge gap of order gA 
the x or radial direction and of order A in the y-direction 
or angular direction. 

To proceed, we must first diagonalize h in Eqn.(130) 
to get new single particle levels, labeled a, b etc. Their 
spacing will also be order A. It is the width of these 
states that we want to obtain, using the self-energy dia- 
gram of Fig. (14). We need couple only a y since a x is 
very heavy. The coupling matrix will change from M x 
to U^M X U, where U relates the a, (3 basis to the a, b ba- 
sis. We will assume that that the rotated matrices have 
the same statistical properties as before rotation, up to 
factors of order unity. The exchanged a y propagator will 
have a mass of order A as can be seen from Eq.(92). 
This results from the fact that motion in the Mexican 
Hat (a change of a y ) has an energy cost of only gA, as 
opposed to a cost of g 2 A in the weak-coupling regime. 
In other words, the Goldstone mode of the self-averaged 
theory acquires a small mass due to sample-specific dis- 
order. The imaginary part of the a y propagator will be 
of the Landau-damped form and we will obtain a T(uj) 
which is essentially the same as Eq.(129). Thus in en- 
tire the strong-coupling regime, the quasi-particle width 
is universal (up to factors of order unity) and of the same 
form as the quantum critical region. This is the meso- 
scopic analog of the Non-Fermi-liquid behavior found in 
the bulk Pomeranchuk phase 44 . 

It is very instructive to analyze Eqn.(128) from differ- 
ent points of view and in different limits, since it describes 
how a very specific and measurable quantity, the quasi- 
particle width, is affected in a real problem where r, g 
and A are at play. 

• (i) First, since in a finite system all many-body 
states are discrete one should compare the typical 
spacing of many-body states with another energy 
scale (such as T, or the level width due to coupling 
to the leads, or the finite energy resolution of the 
experiment, see Silvestrov 38 ) to decide whether the 
width of Eq. (128) is apparent as a continuous 
spectral function or as a set of sharp peaks with an 
envelope given by Eq. (128)^ . 

• (ii) The crossover between the weak-coupling and 



K To compute the decay width at nonzero T we should use 
the formula 

JroE ret ,«,(w) = -g 2 A 2 J2 al3 \M a p\ 2 ± Y, I dw'/m(G ret ,a(w')) 
6(u/ - uj)[N b (uj' -u) + N F (cu')]Im(D re t(cu' - w)) + 



6(w - w')[l + N B (u- w') - N F {u')]Im{Dret{uJ - w'))J (131) 

where Nf is the Fermi function and Nb(uj) = {eP" — 1) _1 is 
the Bose function. 



quantum-critical regimes can be seen in Eq. (128). 
Note that the decay width could have depended on 
the dimensionless parameters r and g in an arbi- 
trary way, but in fact does so only through the en- 
ergy scale £ = 2rgA/ir. This behavior is character- 
istic of systems near a critical point; the right-hand 
side of Eq. (128) is a scaling function, and the com- 
bination £ is a scaling variable near the quantum 
phase transition. Since gA = Et, we see that as 
one approaches the critical point r = 0, the thresh- 
old energy for quasiparticle broadening e* decreases 
as tEt- We emphasize that it never reaches zero; 
at r — 0, it assumes a value A, the infrared cutoff. 
Even if one is in the weak-coupling regime, one can 
probe the physics of the quantum-critical regime 
by looking at the quasiparticle width at the energy 
scale e ~ C — rgA. In past work 37 ' 38 , this behavior 
was expected at ~ y/gA. Thus when r falls below 
l/yfg~, the threshold for such behavior (ill-defined 
quasiparticle) also drops linearly with r to a scale 
as low as A. In units of E T , it drops from E T /y/g 
to E T /g. 

Since our system has three distinct regimes of be- 
havior, it is fair ask how one is to navigate around 
the phase diagram. One option is to alter the par- 
ticle density and hence r s . The above discussion 
gives us another option: To vary oj or e. By increas- 
ing this variable (say by conducting a finite-bias 
study of transport 60 ) we can pass from the weak- 
coupling regime to the quantum-critical regime. By 
contrast, in the strong-coupling regime, the width 
is always as in the quantum-critical regime. The 
reason is that the would-be Goldstone mode a y has 
a very small mass (of order A) entirely due to the 
sample-specific disorder. This is evident from the 
hamiltonian of Eqn. (92) describing fluctuations. 

• (iii) Let us contrast the underlying physics of the 
weak-coupling and quantum-critical regimes. In 
the weak coupling limit, ( = 2rgA/-K is very large 
and we regain 
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a result we could get directly by expanding the <x- 
propagator as 
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and carrying through the ^'-integration in Eqn. 
(122). It is clear that taking higher-order terms 
in the expansion in powers of u m lets the incom- 
ing particle decay into final states with higher and 
higher numbers of particles/holes. From our pre- 
vious large- N analysis we know that the radius of 
convergence of this expansion is since there 

is an instability of the system for \u m \ > |u£J. 
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At the critical point the incoming particle couples 
with equal ease to many-body states with arbitrar- 
ily many particles and holes. 

• (iv) The function in Eq. (129) has no suppression 
with inverse powers of g, in contrast to the weak- 
coupling expression s 2 /g 2 . This is a result of the 
<r gap dropping from order gA to A as explained 
in the points made above. On the other hand MS cX 
function of e it grows quite slowly. The logarithmic 
growth of width with energy justifies the replace- 
ment of the fcrmion spectral function by the delta 
function (Eq. (123)) since the function it multiplies 
in the integral is very broad. 

• (v) Recall that in writing down Eq. (127) we as- 
sumed that the particle-hole states formed a con- 
tinuum even at low energies. Armed with the de- 
cay rate of Eq. (129) we can justify this self- 
consistently. The self-consistent D ret is given by 
the inverse of Eq. (121) with the appropriate self- 
energies put in for e a and sp. Thus, each sharp 
pole in -D^, formerly at e a — ep will now be 
smeared out over an energy range of order T a + Tp. 
If this width is greater than the typical spacing of 
two-body states (which is 2A 3 /e 2 ) the continuum 
approximation is valid. It is easy to verify that this 
holds for energies as low as a few A (modulo the 
caveat in point (i)). 

• (vi) In the quantum-critical and strong-coupling 
regimes we have argued that single-particle states 
are strongly connected to many-body states at 
low energies. One signature of this would be the 
statistics of the many-body states. When a sin- 
gle Slater determinant forms a good description of 
the ground state, one expects the single-particle 
levels to obey Wigner-Dyson statistics, while the 
low-energy many-body states are nearly uncorre- 
cted and obey Poisson statistics. However, in the 
quantum-critical and strong-coupling regimes we 
expect the many-body states at low energies (a 
few A but independent of g) to be strongly cou- 
pled and therefore obey Wigner-Dyson statistics. 
In fact, an exact diagonalization study 63 is con- 
sistent with this picture. Berkovits and Avishai 63 
find that at small coupling strength the many-body 
states display Poisson statistics, while at interme- 
diate coupling strength they obey Wigner-Dyson 
statistics. They also find a reentrance of Poisson 
statistics at very high couplings, possibly a signa- 
ture of a Wigner glass regime 68 ' 69 . 

The decay width of the quasiparticles at nonzero en- 
ergy above the Fermi energy should be visible in non- 
linear conductance experiments 60 (conductance at finite 
bias). 

Finite-T effects in the weak-coupling regime are worth 
a little more elaboration. Assume that one is in the 
weak-coupling regime, whose physics is dominated by the 



Universal Hamiltonian 20 . However, there might be ob- 
servable effects associated with quantum and/or thermal 
fluctuations of er on the peak-height distribution. The 
experiment 50 shows a narrower distribution than that 
expected for nonintcracting electrons. Many ideas have 
been offered to account for this discrepancy, including 
dephasing 64 , the effect of completely random two-body 
interactions 65 , spin-orbit interactions 66 , and finite tem- 
perature and exchange 23 ' 24 . While all these may be con- 
tributing factors, the analysis of Usaj and Baranger 23 
and Alhassid et aP 4 including the effects of exchange 
in the Universal Hamiltonian 20 sense, and the effect of 
nonzero T, seems to be the most comprehensive. How- 
ever, despite the fact that a larger value of the exchange 
constant is used 23 than predicted by RPA (which makes 
the distributions narrower) the experimental distribution 
is even narrower than the prediction, and the discrepancy 
increases with increasing T. 

In our approach, as T increases, higher energy states of 
both the fermions and er are populated, which will modify 
the imaginary part of the fermion Green's function (Eq. 
(131)). This in turn can serve as the input to the formula 
derived by Meir and Wingreen 67 for the current through 
a fully interacting region, which in the simplest case of 
equal couplings to the right and left leads reduces to 

I = ~ J de(N F . L (e) - N FtR (e))Im(Tr[TG}) (134) 

where Np are the Fermi functions of the two leads, T — 
TlTr/(Tl + Tr) is the level width matrix and G is 
the matrix single-particle Green's function (including all 
many-body contributions) . 

It is qualitatively clear that as one approaches the 
quantum critical regime by increasing T or decreasing 
£, the single-particle fermionic states get more and more 
strongly coupled to many-body states. This will lead to 
incoherent averaging of the conductance and thus a nar- 
rowing of the peak height distribution beyond that im- 
plied by the Universal Hamiltonian treatments 23,24 . A 
quantitative calculation is beyond the scope of this pa- 
per. 



C. The Order Parameter and its Fluctuations 

We are now in a position to examine the size of the 
order parameter and its fluctuations in the three different 
regimes, and to delineate the crossovers between them. 
Once again, even and odd m turn out to be significantly 
different, and we will treat them in turn. 

In the case of even m a term linear in er is allowed, 
and the effective potential is generically of the form 

Vefficr) = -go Acos( X - Xo) + ^ 2 Ar<r 2 + ^.g 2 A(<r 2 ) 2 + . . . 

(135) 

where the first term arises from the realization-specific 
disorder and the other terms are self-averaging. The 
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cos(x — Xo) is intended to show the explicit breaking 
of the rotational invariance in er-space by the sample- 
specific disorder. By keeping only the terms kept above, 
and minimizing in the magnitude and angle of er, it is 
easy to verify that near the critical point (r < 1) the 
average value of the order parameter is the solution to a 
cubic and has the form 



and thus 
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where f(x) is a scaling function which becomes a constant 
as x — > and behaves as \/\x\ for large negative x. Thus 
the order paramter behaves as </ -1 / 3 in the quantum- 
critical regime and crosses over to a \f\r~\ behavior deep 
in the strong-coupling regime. From Eq. (136) it is 
clear that the crossover between the weak-coupling and 
quantum-critical regimes happens for |r|g 2//3 w las far as 
the order parameter is concerned. This applies equally to 
the quantum-critical to strong-coupling crossover. Note 
that this is a different crossover than that found in quasi- 
particle broadening. 

For odd m, only even functions of er are allowed in 
V e ff, even among the realization-specific terms. The ef- 
fective potential is generically 

V eff {a) = - 5 Aa 2 cos2( X - xo) + \g 2 Ara 2 + \<?*{o*? 

(137) 

For simplicity, we will assume that either a small nonzero 
external flux selects one of the two degenerate minima 
in the following. Clearly, the crossover to the quantum- 
critical regime happens when the self-averaging quadratic 
term is comparable to the realization-specific term, lead- 
ing to the criterion rg w 1. In the quantum-critical 
regime the order parameter size is of order 1 / y/g, while 

deep in the strong-coupling regime it is behaves as \/\r\- 
Let us now look at the fluctuations of the order param- 
eter. In previous work 10 two of us had wrongly stated 
that quasiparticle broadening was connected to the de- 
localization of a. As has been discussed in the previous 
subsection at some length, quasiparticle broadening in 
the quantum-critical and strong-coupling regimes is a di- 
rect result of the small gap in the er propagator (of order 
A). However, the fluctuations of er remain bounded and 
nearly the same in all three regimes. To see this, let us 
calculate 
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where D is the Euclidean propagator of er. Here wc have 
cut off the frequency integration at E T = gA, because 
retaining higher frequencies would be inconsistent with 
our starting point, in which all states of energy higher 
than Et have been integrated out. From Eq. (124) we 
obtain 
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In the weak-coupling regime £ ~ rg, and the size of the 
er fluctuations is 1/ y/g. This increases a little in the 
quantum-critical and strong-coupling regimes (where £ ~ 
A) to y/log g/g. To reiterate, Fock-space derealization 
has to do with the gaplessness of er, and has no relation 
to its fluctuations. 



D. Summary of Signatures and Connection to 
Experiments 

Since this section has been rather long and involved, it 
is appropriate to collect in this subsection the signatures 
of the existence of the collective variable in the various 
regimes. We remind the reader that we have restricted 
attention to charge-channel instabilities for spinful elec- 
trons, in which the spin index is a passive spectator. This 
restriction is not unreasonable in the light of the values of 
the Landau parameters 42 for r s = 5 (Eq. (21)). We em- 
phasize that only the signatures pertaining to Coulomb 
-Blockade are discussed in detail in this paper. Other 
equally striking signatures in the persistent current will 
be explored exhaustively in Part II 11 . 

• Strong-coupling regime for all m: er is localized in 
a single minimum for m even, or a pair of minima 
related by time-reversal for m odd, with the size 
of |er| being of order unity. The addition spectrum 
and therefore the conductance peak-spacing distri- 
bution are broader and more symmetric than the 
GOE level spacing distribution. This results from 
rare large shifts of er as particles are added, which 
also leads to a correlation between small peak- 
spacing and exponentially small peak height. The 
m even case has a diamagnetic persistent current 12 
in response to an external flux, while the m odd 
case has a paramagnetic persistent current. Low- 
energy quasiparticles are broadened throughout the 
strong-coupling regime, as will be elaborated below. 

• Strong-coupling regime for m odd, isolated dot: At 
T = 0, the system is in a symmetric superposition 
of two exactly degenerate minima of the effective 
potential. However, the crossover flux is exponen- 



tially small in g, of the order of < 
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The breaking of time-reversal by this tiny flux can 
be experimentally monitored in the peak-height 
distribution 39 ^ 40 , which corresponds to the GUE for 
> Across- In the flux regime <p /y/g > 4> > Across 
the peak position varies linearly with <f>, in contrast 
to the noninteracting case in which it would vary 
quadratically. 

• Strong-coupling regime for m odd, dot coupled 
to leads: The dynamics of er becomes dissipative 
at arbitrarily low energies, and our problem can 
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be mapped on to the Caldeira-Leggett problem 3 . 
From its solution 35,36 , we conclude that at suf- 
ficiently high dissipation, corresponding to level 
widths (due to the dot- lead coupling) of T ~ A/y/g, 
the system can be driven to spontaneously break 
time-reversal symmetry. This transition can once 
again be experimentally accessed via the peak- 
height distribution 39,40 , and is isomorphic to a fer- 
romagnetic to antiferromagnetic Kondo transition. 

• Quasiparticle width: Throughout the quantum- 
critical and strong-coupling regimes er becomes 
"gapless" in the scaling limit (g — > 00, A — > 
with Et hxed). A single-particle description is 
no longer valid even at low energies (of order a 
few A), and low-energy quasiparticles are "Fock- 
space delocalized" . This prediction is in contrast to 
previous work 37,38 which suggested a broad quasi- 
particle only at energies of order ~ A^/g. Physi- 
cally this corresponds to strong mixing between the 
low-energy single-particle excitations and collective 
overdamped excitations of er. In the quantum- 
critical regime both vector components of er are 
gapless, while in the strong-coupling regime, the 
fluctuations of er in the nearly degenerate Mexican 
Hat potential (the would-be Goldstone mode) are 
nearly gapless. 

• Weak-coupling regime: Despite the fact that it 
is dominated by Universal Hamiltonian physics, 
quantum and/or thermal fluctuations of er might 
have a measurable effect on the width of the peak- 
height distribution. Such effects are predicted to 
get stronger as one goes towards stronger coupling 
(by increasing r s ). The quantum-critical regime 
can be accessed by measurements at hnite bias, 
nonzero frequency, or nonzero T. 

Let us now briefly analyze the relevant experiments. 
The Sivan et al 5 and Patel et aP experiments are done 
on gated GaAs 2DEG samples, for which r s ss 1 — 1.2 in 
the bulk. Using the area of the sample and the fact that 
these dots are in the ballistic limit we can find both A 
and Et, and thence find j«7- 14. Sivan et af find that 
the CB peak spacing is about 5 times broader than that 
predicted by the Universal Hamiltonian, which describes 
the weak coupling region of our phase diagram. However, 
Patel et al 5 find it to be in accord with this prediction, 
after accounting for "experimental noise" which is de- 
termined by measuring the magnetic field asymmetry of 
the CB spacings (which indicates motion of the dopant 
atoms, or some other scrambling of the single-particle 
potential). Thus, the Patel et al data seem to lie in the 
weak-coupling region. The Sivan et al data also corre- 
spond to the weak-coupling region. Presumably, if the 
subtraction of experimental noise were to be carried out 
this data would also match the Universal Hamiltonian 
predictions. 

The experiments by Simmcl et al and Abusch-Magder 
et af are performed on Si quantum dots with r s w 2.2, 



and g s=s 18. (For r s > 2 local charge density correla- 
tions develop in the dot 68 , similar to the classical limit 69 . 
While a Fermi surface distortion is not a charge density 
wave, it enhances the susceptibility for one, and could 
thus be a precursor). Two signatures of the critical fan 
are found in this experiment 6 . The CB peak-spacing dis- 
tribution is found to be 7-8 times wider than expected 
from the Universal Hamiltonian (assuming spin degener- 
acy). Also, the width of the CB peaks does not vanish 6 
as T — > 0. This is just what is expected for a system 
was located in the critical fan: The ground and low-lying 
excited states are "Fock-space delocalized", and single- 
particle states are broad even at low energy. However, it 
is not clear whether the width of the peaks is intrinsic or 
dominated by charge rearrangement noise 70 . 

Finally, as has been discussed in the previous sub- 
section, measurements of the width of the peak-height 
distribution 50 and its evolution with T arc not com- 
pletely understood from the Universal Hamiltonian 
standpoint 23 ' 24 . Our theory predicts a further narrow- 
ing of the peak-height distribution. 

In conclusion, current experiments seem to be in 
the weak-coupling regime, with possible excursions into 
the quantum-critical crossover at nonzero temperatures. 
From the values of the Landau parameters for the clean 
2DEG 42 one can guess that systems with r s substantially 
larger than 5 will be needed to see the strong-coupling 
regime, and that charge-channel instabilities are likely 
to occur first, although occasionally even systems with 
r s = 5 may happen to be in this regime (see the end 
of Appendix A). In constructing more strongly interact- 
ing sample, one must bear in mind that since g ~ kpL, 
reducing the density at constant dot size L will reduce 
g. What one would ideally like to do is to keep g fixed, 
which means preparing larger dots while reducing elec- 
tron density. 

VI. CONCLUSIONS AND OPEN QUESTIONS 

Starting with the problem of interacting two- 
dimensional electrons in a ballistic mesoscopic structure 
with chaotic boundary scattering, we found that Landau 
Fermi-liquid interactions (parameterized by dimension- 
less Landau parameters u m in the m th angular momen- 
tum channel) were the most natural starting point at the 
Thouless scale Et- The only condition for this to hold 
is that Ef/Et be large. For a ballistic/chaotic quantum 
dot this ratio is of the same order as the dimensionless 
conductance g = Et/A, where A is the mean level spac- 
ing. Thus the largeness of g is a sufficient condition for 
Fermi liquid interactions to be applicable in the Thou- 
less band. This is a new regime in the broad problem 
of the interplay of interactions and disorder. The special 
features of this regime which make it controllable are the 
applicability of Random Matrix Theory, and the pres- 
ence of the small parameter 1/g, which conspire to make 
a large- N (N — g) approximation feasible. The solution 
is nonperturbative in both interaction strength and dis- 
order: The disorder is the strongest it can possibly be, 
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scattering momentum states chaotically into each other, 
while the interactions have a dimensionless strength of 
order unity. The secret to the success of the nonper- 
turbative solution is the additional small parameter 1/g 
which controls the size of fluctuations in the large- N ap- 
proximation. Note that in contrast to many applications 
in which the large- N approach is used, in this case it is 
clearly justified by the largeness of g (5-20 in experimen- 
tal samples). 

Restricting our attention to charge-channel instabili- 
ties, we found a disordered version of the Pomeranchuk 
shape transition in every Landau Fermi-liquid channel, 
confirming the one-loop RG result found earlier by two 
of us 9 , and putting it on a rigorous footing as being ex- 
act in the large-g limit 10 . The corresponding Landau pa- 
rameter had to be u m < u* n = —1/2 In 2 = —0.7213, as 
compared to the bulk instability at u m < — 1. Thus there 
is a window in which the clean bulk system is stable, but 
the mesoscopic ballistic/chaotic system is not, which in- 
dicates that both interactions and disorder are crucial 
to the existence of this regime. More specifically, we 
saw that interactions produce a Mexican Hat potential, 
and the realization-specific disorder allowed symmetry- 
breaking even in the presence of finitc-g quantum fluctu- 
ations. This can be understood as explicit mass genera- 
tion for the would-be Goldstone mode by sample-specific 
disorder. 

The transition at g — oo is replaced by a crossover 
for any finite g. In accord with generally known features 
of quantum phase transitions 31 the coupling constant- 
ly plane is divided into a weak-coupling regime, a 
strong-coupling regime, and a fan-shaped quantum crit- 
ical regime. However, our theory displays additional 
nontrivial features not found in clean models, such as 
the persistence of symmetry-breaking even in the zero- 
dimensional limit, due to sample-specific disorder. In 
particular, in our theory the shape of the crossover 
regime depends on the property being considered: For 
example, the quasiparticle width at low energies has 
only a single crossover between the weak-coupling and 
quantum-critical regimes and remains broad throughout 
the strong-coupling regime. On the other hand, the or- 
der parameter displays an additional crossover from the 
quantum-critical to the strong-coupling regime, with the 
shape of the crossover depending on whether the angular 
momentum channel m in which the instability occurs is 
even or odd. It must be reiterated that these proper- 
ties can be explicitly calculated in a controllable large- TV 
approximation in all the regimes thanks to the small pa- 
rameter 1/g. 

Note that we make no attempt to keep all 1 j g contri- 
butions in the effective action S e f / , instead keeping only 
the lowest nontrivial terms of each type. As an exam- 
ple, we do not keep 1/g corrections to the self-averaging 
part of the effective action, but we do keep the 1 /g terms 
which break the rotational invariance in er-space of the 
self-averaged part, since they have a profound effect on 
the physics of the strong-coupling regime. Similarly, we 
keep the "kinetic" terms in S e f / which are down by 1/g 
compared to the self-averaging static part of 5 e //, since 



these are the lowest order nonzero terms. 

The strong-coupling regime is where the most striking 
properties of the collective Fermi surface distortion field 
cr are found. Here a bifurcation occurs with some proper- 
ties being common to both even and odd m, while others 
depend on m modulo 2. The features shared by all m in 
the strong-coupling regime include: (i) A mean value of 
the order parameter of order unity, (ii) Rare large shifts 
of the order parameter upon adding a particle, leading 
to a broadened peak-spacing distribution, (iii) The same 
phenomenon as (ii) leading to a correlation between small 
peak-spacing and exponentially small (~ exp— ghig) 
conductance due to the near-orthogonality 61 ' 62 of the 
states at different (cr). (iv) Quasiparticlcs become broad 
at energies of a few A throughout the strong-coupling 
regime. A special property of the m even strong-coupling 
regime is a diamagnetic persistent current 12 , to be inves- 
tigated in greater detail in Part II 11 of this paper. 

The strong-coupling regime of a system with an in- 
stability in an odd m channel shows even more striking 
effects, which depend sensitively on how strongly the dot 
is coupled to the leads. The physics can be mapped on to 
the Caldeira-Leggett model 34 , and is driven by the com- 
bination of an exact two-fold degeneracy in the effective 
potential (a consequence of the time-reversal invariance 
of the underlying Hamiltonian) and the low-energy dis- 
sipation induced by coupling to the leads. At weak dot- 
lead coupling, the system is "almost" in a time-reversal- 
(T-) broken state. The crossover flux need to tip the 
system into a completely T-broken regime is exponen- 
tially small (Across — 4>oe~ 9 /g), as contrasted with a 
crossover flux of 0o / \fg needed for a noninteracting sys- 
tem. When the dot-lead coupling increases to a value 
such that the broadening of single-particle states due to 
the leads reaches T ~ A/^Jg, a spontaneous breaking of 
time-reversal takes place 35,36 . This is a true quantum 
phase transition in the universality class of the antifcr- 
romagnetic to ferromagnetic Kondo model 55 , and can be 
acheived by a finite system by coupling coherently to an 
infinite reservoir. In the T-broken phase the system also 
displays a spontaneous paramagnetic persistent current 12 
of order Et/4>q, to be explored in greater detail in Part 
II 11 . 

Coming to the quantum-critical regime, the main fea- 
ture is the "gaplessness" of the collective variable cr, 
by which we mean that its propagator has the Landau- 
damped form — 1/ioj for A <C w <C Et- In this regime 
low-energy single-particle excitations are strongly cou- 
pled to many-body excitations with arbitrary numbers 
of particle- hole pairs via the collective variable cr. This 
has many consequences: (i) Recall that the threshold 
for a continuous and broad spectral function for a parti- 
cle was found to be e* ~ Ay/g/hig in earlier work 37 ' 38 
which ignored correlations between interaction matrix el- 
ements, and therefore missed the quantum phase tran- 
sitions found here. In our approach we find that this 
threshold e* decreases to a few A independent of g 
in the quantum-critical regime. This could potentially 
be seen in nonlinear conductance (conductance at fi- 
nite bias) experiments 60 , (ii) Since they are strongly 
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mixed by interaction effects, the statistics of the low- 
energy many-body states should change in the quantum- 
critical regime from Poisson towards Wigner-Dyson, as 
has been seen in recent numerics 63 (iii) Finally, and per- 
haps most relevant for current experiments, one should 
be able to access the quantum-critical regime from the 
weak-coupling regime by finite-bias, finite-frequency, or 
finite-temperature measurements. In particular, there 
should be a narrowing of the peak-height distribution be- 
yond Universal Hamitonian effects 23 ' 24 at nonzero tem- 
perature. 

As has been mentioned above, the broadening of 
low-energy quasiparticles persists in the strong-coupling 
regime due to the near-gaplessness of the would-be Gold- 
stone mode of the Mexican Hat. This seems to be 
the mesoscopic analog of the non-Fcrmi-Liquid behavior 
found throughout the bulk Pomeranchuk phase in clean 
systems 44 . 

Finally, as regards the weak-coupling regime, the col- 
lective variable <x has a huge gap here of order rg (r = 
Ju~\ ~ TtFT 1S ^ e distance from the critical point). Thus 
the weak-coupling regime is controlled by the Universal 
Hamiltonian 20 of Eq. (3), which is the low-energy fixed 
point in this regime. 

The picture that emerges is that just as Landau's 
Fermi liquid theory is unstable to the bulk Pomeranchuk 
instability in any Fermi liquid channel, the Universal 
Hamiltonian is also unstable to the mesoscopic Pomer- 
anchuk instability in any channel, with the instabil- 
ity occuring more readily in mesoscopic ballistic/chaotic 
systems than in the bulk. The mesoscopic Pomer- 
anchuk regimes represent physics nonperturbatively dif- 
ferent from Universal Hamiltonian physics. These strong- 
coupling regimes, and the crossover quantum-critical 
regimes, can be accessed by changing the electron den- 
sity, or temperature, or the probe frequency. The con- 
silience of Random Matrix Theory and large-g enables 
us to carry out controlled calculations and make predic- 
tions in all three regimes. Our analysis provides a new 
framework for analysing future experiments on strongly 
interacting mesoscopic systems. 

Let us now turn to open questions and possible fu- 
ture directions. Apart from a more precise treatment 
of finite- T effects 22-24 in our theory, the most press- 
ing open question is the treatment of spin-channel in- 
stabilities, which necessarily involves the inclusion of 
the Universal Hamiltonian exchange coupling 20 J. In 
other words, what is the interplay between the meso- 
scopic Stoncr transition 20 and the mesoscopic Pomer- 
anchuk transition in a spin channel? One can imagine a 
cooperative phenomenon whereby the exchange coupling 
favors symmetry-breaking in a spin Fermi-liquid channel. 
This might reduce the magnitude of the critical coupling 
considerably, and hence make the transition observ- 
able at smaller r s . Exotic types of ordering intertwining 
spin and momentum space indices might also occur 43 . 

A second effect which deserves investigation is the fol- 
lowing: We have assumed that no renormalization of the 
Landau parameters happens at energy scales higher than 
Et, whereas they do renormalize below Et in the man- 



ner described in Section III. Clearly, there must be some 
renormalization of the Landau parameters above Et, 
and since the sign of the flow cannot abruptly change, 
this effect also tends to reduce the critical coupling. To 
rephrase the question, are the Landau parameters of the 
chaotic/ballistic cavity the same as those in the bulk 
2DEG of the same electron density? Perhaps the best 
way of resolving this issue is to carry out the analog of 
the bulk quantum Monte Carlo calculations 42 directly in 
a chaotic/ballistic cavity 71 . 

A somewhat related question concerns slightly diffu- 
sive quantum dots, where the bulk mean free path I is 
smaller than (but not negligible compared to) the size 
of the dot L. In this case there is an energy regime be- 
tween Tivf/L and the Thouless energy Et = fiD/L 2 . In 
this intermediate energy regime there will be a crossover 
between bulk diffusive dynamics and Random Matrix dy- 
namics. In the extreme diffusive limit one must recover 
the Finkclshtcin scaling 45 ' 46 . 

One can also ask what happens as one increases the 
interaction strength (or r s ) to very large values. Is 
there physics beyond the strong-coupling regimes we have 
found? Certainly for the Coulomb interaction we know 
that one must eventually cross over into a Wigner Glass 
(or Coulomb Glass) regime 68,69 . The flattening of parts 
of the Fermi surface in the mesoscopic Pomeranchuk 
regimes increases the susceptibility for charge density 
wave instabilities. This indicates that the crossover to the 
extremely strong-coupling Wigner Glass regime might 
happen more readily in mesoscopic ballistic/chaotic sam- 
ples than in bulk disordered samples. An investigation 
of this crossover would complete the picture of strong- 
coupling physics in mesoscopic ballistic/chaotic samples. 

At a more philosophical level, the identification of a 
new, controllable, regime of disorder and interactions 
opens up many possibilities. As an example, one could 
apply our RG techniques with RMT correlations to the 
Kondo problem 55 , and investigate the effect of Landau 
Fermi-liquid interactions between conduction electrons 
residing in a ballistic/chaotic cavity 72 on the Kondo ef- 
fect. Another example to which our techniques might 
be applied is the effect of Fermi-liquid interactions on 
granular, disordered, gapless superconductors exhibiting 
a novel metal-insulator transition 73 , since the eigenvalue 
and eigenvector statistics of their grains are expected to 
be governed by one of four newly discovered RMT uni- 
versality classes 74 . 

To conclude, ballistic/chaotic quantum dots seem to 
provide us with a unique theoretical and experimental 
playground with exquisite control on both fronts. On 
the theoretical front, RMT and the 1/g expansion allow 
us to tackle disorder and interactions simultaneously and 
nonperturbatively. On the experimental front dots give 
us the best of both worlds: A phase transition at g = oo 
which we can "see" even at finite g and access to all the 
control parameters a mesoscopic system allows (varying 
r s , external flux, bias voltage, and coupling to leads). 
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VII. APPENDIX A 

In this appendix we will show that the RG equation 
and the quadratic part of the effective action are self- 
averaging. Let us first focus on the RG, and look at the 
part of Eq. (40) which is to be summed over 



R= E 3qfcT£«(0k-0p)«OV-0k') 



-A 1 " 1 pp' 

x C(p)C(p')<Mp)0m(p') 



(141) 



As has been stated before, in RMT the wavefunction av- 
erages do not depend on the energy separations of the 
states, so the wavefunction average can be carried out 
separately. This average in the GOE for fi ^ v and 
generic momentum labels is Eq. (41), reproduced below 
for convenience 



< ^(Pl)</C(P2)<MP 3 )^(P4) > = 
<5l4^23 <$13<524 + <*>!, -2^3,-4 



(142) 



The 1/g 2 term is the "naive Wick contraction" of 
leading-order RMT, but the 1/g 3 term is necessary to 
maintain orthogonality between /i and v. The final term 
would be missing in the GUE. Substituting the correct 
momentum labels for the particlc-holc diagram we see 
that the wavefunction average is 



3pp' 



1 + 6, 



p ,-p' 



(143) 



Let us first carry out the ensemble average of R. Using 
Eq. (143) in Eq. (141) we obtain a convolution of the 
two Fermi liquid functions 

$>(f?k - p )u(9 p -6*)= g{u 2 + ~Y, u2 m ^sm(6 - 9')) 



m—1 



(144) 



where we have reverted to the notation 9 = 9k, 9' = 6*k' • 
Note that we have used 



d9p 

2tt 



(145) 



which is valid for large g. In the second term of Eq. 
(143), the <5p,- P ' turns out to be subleading, while the 
other allows independent sums over p, p'. This means 
that only uq contributes due to this term, which produces 



5>(*k - 0p)u(0 p , -*k') =9 
pp' 



2 u 2 (146) 



Next we need to do the energy sum, which we replace 
by an integral 



51 TXT 



£„=-A 



ds 1 In 2 



A +\e v \ J AA + e A 



(147) 



Putting all the pieces together, we obtain the ensemble 
average of R to be 



« i?> = 



In 2 



a m=l 



U m COS TO 



9') 



(148) 



Now consider the fluctuation of the quantity R of Eq. 
(141), that is, R— < Each term in R is of or- 

der 1/g 3 , where one factor of 1/g comes from the energy 
denominator, while each wavefunction contributes 1/ ^fg. 
There are g 3 such terms coming from the three summa- 
tions. Since the average has been taken out, each contri- 
bution is completely random and comes with a random 
sign. From a one-dimensional random walk argument, 
with each step being of size 1/g 3 , and g 3 steps, the net 
answer must be of order 



R- <i?> 



1 

^3/2 



« < R > 



(149) 



Thus we see that the average of R indeed dominates 
its fluctuations over the ensemble, and R is therefore self- 
averaging. 

The situation is simpler for the quadratic term of the 
effective action, Eq. (79) (call it S^). 



N F (0) - N F (a) 



^C(k)^(k)^(k')i(k') 



a/3 £a ~ ^ k,k' 

(a 1 cos Tti9 + cr 2 sin mO) (01 cosm9' + a 2 smm9') (150) 

In effect, S2 is the same as R, except for the fact that 
(i) only a single Fermi liquid channel contributes, and (ii) 
there is an extra sum over /i, which was not present in R. 
The self-averaging of S2 follows the same lines as above. 
Because of the extra sum (and the factor of g 2 in front) 
the ensemble average <C S2 3> is of order g 2 . Considering 
the same random- walk argument for £2— -C S2 S> we 
have four sums (g 4 steps in the random walk) with each 
step being of order 1/g 3 , and an additional factor of g 2 
in front, leading to 



S 2 - < > ~ g 



~ g « < S 2 > (151) 



Thus the fluctuations of S2 are 1/g down from its self- 
average. 

There is one subtlety with the effective action which 
is not present in the RG equation. The denominators 
of the RG equation can never be smaller than A, while 
those of 5*2 go all the way to the smallest gap near the 
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Fermi energy. In the GOE the ensemble average <§C 5| ^> 
formally diverges due to the linear behavior of the level 
spacing distribution for small spacings. This implies that 
while the distribution for S 2 is sharply peaked around its 
average, it also has a long tail, which implies that the 
critical coupling has a long tail as well. Calling the 
change in the critical coupling 5u* , and taking only the 
smallest level spacing into account, we find the tail to be 



Ptaii(Su*) w l/g 2 (5u*f 



(152) 



This has the consequence that even systems nominally 
in the weak-coupling regime may occasionally be in the 
strong-coupling regime, with a probability that goes as 
I / g 2 (5u*) 2 . Considering a system at r s = 5, with $ 2 = 
—0.25, we have 5u* = 3.279. Assuming j = 5we have a 
probability of 0.0037 of the system being in the strong- 
coupling phase. 



VIII. APPENDIX B 

Consider a generic theory with a fermion and a boson 
(ip and er in our problem), in which we integrate out 
the fermions to get an effective action S e ff for er. It is 
clear that this action gives all the answers for questions 
involving er alone. However, if at this stage we want to 
ask a question involving the fermions, which have been 
integrated out, how should one proceed? This is not an 
unprecedented situation, but we discuss it nonetheless for 
completeness. 

Consider the following schematic partition function 
in which the fermions have been coupled to fermionic 
sources J, J to enable a computation of their correlators: 

Z(J, J) = J dadipdipe s (153) 

S(a, ip, i>) = Jtp + J^P + i>(d + Ma)ip - a 2 /2 (154) 

where M is a coupling matrix, and d denotes the fermion 
kinetic energy. (In our case (d = iu) — £ a ))- 
Upon integrating the fermions out we get 



s eff(<r) = J 



d + Ma 



J + Tr ln(<9 + Ma) - a 2 /2. (155) 



To obtain the fermion propagator G, we take the J 
and J derivatives of In Z, set J = J = 0, and obtain 

G=(Z(0,0))" 1 J rf(7 ^J_- e ^W. (156) 

Thus the full G is just the propagator in a given ex- 
ternal field <t, functionally averaged of over a with its 
effective action in the Boltzmann factor. 

Let S e f f have the following expansion about a saddle- 
point: 

S e ff = So(&o) + -5aD~ 1 (ao)5a + neglected terms 

(157) 



If we ignore the fluctuations we get 



d + Ma a 



(158) 



If we keep the gaussian fluctuations, we get to lowest 
nontrivial order 



G 



1 



d + Ma a 



1 



-M- — 

2 d + Ma 



M 



1 



d + Ma 



(5a da) 



where the average 



(8a 5a) = D(a ) 



(159) 



(160) 



is just the a propagator in the gaussian approximation. 
Clearly, this procedure generates all the Feynman di- 
agrams containing the three-point vertex coupling the 
fermions to fluctuations of er. This series can be re- 
summed in the usual way by constructing a self-energy 
for the fermions, with the lowest nontrivial contribution 
to it being Fig. 15. Thus there is no double-counting 
involved in using Fig. 15 (or Eq. (122)). 

In the weak-coupling regime, a = 0, and Eq. (121) 
gives the inverse a propagator. In the symmetry-broken 
strong-coupling regime, the term d + Mctq corresponds 
to the hamiltonian of Eqn. (130), the treatment of which 
follows the equation. 
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